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ABSTRACT 


A  five  degree  of  freedom  (constant  pitch)  simulation  model  for  a 
Small  Waterplane  Area  Twin  Hull  (SWATH)  ship  has  been  developed.  Complete 
descriptions  of  the  mathematical  models  employed  and  a  description  of  the 
corresponding  digital  simulation  program  are  presented  in  this  report.  Per¬ 
formance  predictions  based  on  the  mathematical  model  are  also  included.  The 
mathematical  models  are  based  primarily  on  experimental  model  test  data 
measured  at  the  rotating  arm  facility  at  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  (DTNSRDC).  The  performance  predictions  gen¬ 
erated  in  this  study  are  appropriate  for  carrying  out  rudder  configuration 
design  evaluations.  The  results  indicate  that  the  spade  rudder  configura¬ 
tion  is  adequate  for  meeting  turn  diameter  specifications. 
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INTRODUCTION 


The  Navy  is  currently  carrying  out  preliminary  design  and  feasi¬ 
bility  studies  on  various  Small  Waterplane  Area  Twin  Hull  (SWATH)  ship 
configurations  in  the  2000  -  20,000  ton  weight  range.  These  ships  have 
desirable  performance  properties  due  to  their  design.  The  SWATH  ship 
basic  configuration  consists  of  two  submerged  hulls  attached  to  two  vertical 
struts.  Propulsion  is  provided  by  propellers  located  at  the  end  of  the  hulls 
Since  the  waterplane  is  small,  surface  waves  have  a  smaller  effect  on  SWATH 
ship  motion  than  on  conventional  craft  which  have  larger  waterplane  areas. 
Also,  SWATH  ships  can  turn  In  a  fairly  small  circle  at  low  speed  by  applying 
a  differential  thrust.  A  sketch  of  the  basic  SWATH  configuration  is  illus¬ 
trated  in  Figure  1. 

In  order  to  carry  out  evaluation  studies,  the  necessary  tools  must 
be  developed.  One  such  tool  is  a  time  domain  simulation  model  suitable  for 
predicting  ship  behavior  during  various  maneuvers.  The  simulation  model  must 
be  sufficiently  flexible  so  that  different  SWATH  configurations  can  be  easily 
evaluated  and  compared.  An  appropriate  simulation  tool  has  been  developed 
during  the  study  reported  here. 

The  general  approach  used  in  developing  the  SWATH  simulation  model 
involves  the  development  of  mathematical  models  based  on  experimental  force 
and  moment  data.  The  experimental  data  is  first  processed  to  develop  a  suffi 
ciently  dense  data  base  by  interpolating  and/or  extrapolating  for  missing 
grid  points.  Also,  symmetry  properties  are  employed  to  generate  a  complete 
data  base  required  to  characterize  the  ship  over  all  orientations  and  speeds. 
The  force  and  moment  values  for  a  given  ship  state  are  obtained  via  multi¬ 
dimensional  interpolation  of  the  data  base.  Theoretical  expressions  are 
used  whenever  experimental  model  data  is  lacking. 
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FIGURE  1  SCHEMATIC  OF  SWATH  6A  CONFIGURATION 


The  resulting  force  and  moment  models  are  used  in  the  differential 
equations  that  describe  the  ship  motion.  These  equations  are  solved  numeri¬ 
cally  using  standard  integration  procedures. 

Various  sets  of  experimental  data  are  used  in  the  present  effort. 
These  data  characterize  the  hull,  rudder,  fin,  and  propulsion  subsystems. 
These  subsystems  are  modeled  separately.  However,  some  interference  effects 
are  taken  into  account.  For  example,  the  interference  of  a  hull  on  its 
propeller  is  accounted  for  by  a  wake  reduction  factor  based  on  experimental 
data.  Other  interference  effects  are  assumed  negligible.  These  assumptions 
are  validated  via  an  analysis  of  the  experimental  data  in  most  cases. 

The  mathematical  model  development  is  presented  in  the  next  section. 
Subsequent  sections  present  results  and  digital  computer  program  documenta¬ 
tion.  Actual  computer  output  and  card  images  of  various  test  runs  are 
presented  in  the  appendices. 
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MATHEMATICAL  MODEL  DEVELOPMENT 


BASIC  EQUATIONS  OF  MOTION 

The  equations  governing  the  motion  of  the  SWATH  are  extensions  of 
Newton's  second  law,  F  =  Ma,  to  the  six  degrees  of  freedom  possible  for  ship 
motion.  The  SWATH  simulation  model  is  based  on  the  standard  equations  of 
motion  for  a  submarine  outlined  in  NSRDC  Report  No.  2510.  Minor  changes  have 
been  made  in  converting  the  submarine  equations  to  a  set  of  equations  appro¬ 
priate  for  a  surface  ship.  The  predominate  coordinate  system  used  in  the 
simulation  model  is  a  body  fixed  system  as  defined  in  reference  #2510.  All 
forces  and  moments  are  defined  in  terms  of  the  body  system.  Unless  otherwise 
specified  all  parameters  are  given  in  the  MKS  system  of  units. 

The  equations  of  motion  presented  in  NSRDC  Report  No.  2510  can  be 
written  in  the  generalized  form, 

^  -l 
S  =  A  1  F 

where  S  is  the  state  vector,  A  is  a  generalized  mass/inertia  matrix  and  F 
is  a  generalized  force  vector. 

The  generalized  mass/ inertia  matrix  appropriate  for  SWATH  is  pre¬ 
sented  in  Table  1.  All  symbols  used  in  Table  1  and  in  subsequent  sections 
of  this  report  are  defined  in  the  notation  section.  The  generalized  force 
vector,  F,  contains  kinematic  terms  plus  contributions  from  the  fin,  hull, 
rudder,  aerodynamic,  buoyancy,  propeller,  and  rate  dependent  forces  and 
moments.  The  components  of  the  generalized  force  vector  are  defined  in 


Table  1  -  GENERALIZED  MASS/INERTIA  MATRIX 


rud  hull  air  Abuoy  prop 

-  w*q)M  +  |(q2  +  r2)x  -  pqy  -  pr2  I 


Table  2  -  GENERALIZED  FORCE  VECTOR 


I 


JM 

m  im 

The  time  rate  of  change  of  the  roll,  pitch,  and  yaw  angles  are 
obtained  by  solving  the  kinematic  equations  presented  in  Report  Mo.  2510 
for  <j>,  0,  and  tp.  The  results  are, 

4>  *  (q  sine  sin«j>  +  r  sine  cos4>  +  p  cose)/cose,  (2) 

0  -  q  cos4>  -  r  sin<}>,  (3) 

4i  *  (q  sin<js  +  r  cos<j>)/cose.  (4) 


Velocity  components  in  the  fixed  or  Earth  coordinate  system  are 
obtained  by  transforming  the  body  x,  y,  and  z  velocity  components  to  the 
fixed  system: 


- 1 

X* 

U 

u 

il 

-il 

V 

w 

The  transformation  matrix  T  is  given  in  Table  3. 

The  above  equations  of  motion  are  suitable  for  a  full  six  degree 

of  freedom  simulation  model.  The  SWATH  ship  considered  in  this  report  has 

four  fins  that  can  be  used  to  provide  a  pitch  moment  such  that  the  pitch  angle 

remains  constant,  i.e.,  0  *  0.0..  Rather  than  solve  the  extremely  complicated 

problem  of  determining  a  fin  angle  of  attack  that  would  result  in  a  value  of 
• 

q  that  would  in  turn  result  in  a  constant  pitch,  the  kinematic  equations  can 
be  used  to  obtain  a  differential  equation  for  q  such  that  the  pitch  remains 
constant.*  When  0  =  0.0,  Eouation  (3)  reduces  to 


*  9  =  0.0  does  not  imply  q  =  0.0  unless  <P  *  0.0  or  r  =  0.0. 
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TABLE  3  BODY  TO  FIXED  COORDINATE  SYSTEM  TRANSMISSION  MATRIX 


A  differential  equation  for  q  is  obtained  by  taking  the  derivative  of  the 
above  equation : 

q  cos4>  -  r  sin<js  =  4>  (q  si n4>  +  r  cos<j>).  (7) 

<j>  can  be  eliminated  in  the  above  equation  by  using  Equation  (2).  The  constant 
pitch  simulation  model  employs  the  above  equation  for  q  rather  than  the  q 
equation  given  in  Report  No.  2510.  Conceptually,  the  above  approach  is 
equivalent  to  including  a  fin  pitch  moment  in  the  pitch  equation  such  that 
9  =  0.0.  In  fact,  if  a  strongly  restoring  pitch  moment  of  the  form 

m  =  ce 

is  included  in  the  differential  equation  for  q  in  six  degrees  of  freedom, 
the  pitch  moment  remains  at  a  value  of  zero  degrees.  Such  a  test  was  conducted 
in  order  to  verify  the  five  degree  of  freedom  equations. 


HULL  HYDRODYNAMIC  FORCE  AND  MOMENT  MODEL 


The  mathematical  models  used  to  represent  hull  forces  and  moments 
are  based  on  experimental  model  data  provided  by  DTNSRDC.  The  data  supplied 
was  obtained  from  experiments  conducted  on  a  1/22.5  size  scale  model  at  the 
rotating  arm  facility  located  at  DTNSRDC.  Six  different  model  configurations, 
as  described  in  Table  4  ,  were  considered.  The  results  of  these  experiments 
were  supplied  in  the  form  of  parametric  curves  for  non-dimensional  values  of 
Y,  K,  and  N  versus  non-dimensional  turn  rate  for  various  drift  angles. 
Measurements  were  not  made  for  isolated  grid  points.  These  missing  values 
were  obtained  by  interpolation  or  extrapolation  based  on  an  inspection  of 
various  data  trends.  The  final  data  bases  used  to  characterize  the  hydro- 
dynamic  hull  forces  and  moments  are  presented  in  Tables  A1  -  A4.  A  four 
variable  table  look-up  or  interpolation  procedure  is  used  in  the  digital 
simulation  model  to  obtain  values  of  Y,  K,  and  N  at  given  values  of  the 
speed,  roll  angle,  drift  angle,  and  non-dimensional  yaw  rate  from  the 
dimensional  form  of  the  data  base.  Interpolations  were  initially  performed 
on  the  non-dimensional  form  of  the  data  base.  However,  due  to  the  fine 
sampling  in  speed  (5  knot  spacings)  linear  interpolation  of  the  dimensional 
form  of  the  data  is  sufficiently  accurate  and  much  more  efficient.  A 
comparison  of  both  results  is  illustrated  in  Figures  2-4.  The  experimental 
data  plotted  in  these  figures  are  hand  faired  using  French  curves.  The  very 
close  agreement  between  the  interpolated  values  and  the  smooth  curves  indi¬ 
cate  the  high  degree  of  accuracy  possible  by  use  of  an  interpolation  pro¬ 
cedure.  Note  that  the  data  bases  contain  rudder  contributions  at  zero  rudder 
angle.  Also,  it  is  assumed  that  the  force  coefficients  can  be  superimposed 
due  to  the  strong  linear  nature  of  the  forces  and  moments. 


TABLE  4  DESCRIPTION  OF  TEST  CASES  A 


NUM-DiMhMSIONAL  Y  (x  10 


niMFNSlONAl 


NON- DIMENSIONAL  K  (x  10 


5.0 


5.0 


DIMENSIONAL  K(x  10 


NON-DIMENSIONAL  N  (x  10' 
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DIMENSIONAL  (x  10M 


81/iuw'dJ  if 


The  range  of  variables  covered  in  the  final  data  base  used  in  the 
mathematical  model  proved  adequate  for  all  simulated  maneuvers  except  for 
turns  carried  out  by  application  of  a  differential  thrust  at  5  knots.  For 
these  cases,  the  non-dimensional  turn  rate  (r1  X  r/u,  u-O)  can  be  as  large  as 
1.5  and  drift  angles  as  large  as  10.0°.  However,  a  linear  extrapolation  of 
the  available  data  should  produce  reasonable  results  since  all  low  speed 
parametric  curves  examined  show  strong  linearity. 

The  drag  alone  the  body  x  direction  is  calculated  from  EHP  curves 
provided  by  DTNSRDC.  For  a  given  value  of  the  EHP,  the  drag  is  obtained  from 
the  formula, 

Orag  =  (EHP  •  326 )/U, 

where  U  is  the  speed  in  knots.  This  mathematical  is  strictly  valid  for 
straight  ahead  flight  at  zero  trim.  Effects  of  ship  motion  can  be  easily 
taken  into  account  when  more  extensive  EHP  data  becomes  available.  The 
effect  of  changes  in  drag  due  to  changes  in  ship  immersion  is  taken  into 
account  by  computing  the  percent  changes  in  frontal  and  wetted  areas  from 
ship  geometrical  factors.  The  resulting  drag  formula  is  given  by, 

R'  =  R( 1.0  +  0.08(Zq  +  Z£)), 

where  R  is  the  nominal  resistance  that  would  be  obtained  for  Z^  =  0.0,  ZQ 
is  a  reference  heave  in  meters,  and  Z^  is  the  heave,  which  is  assumed  posi¬ 
tive  downward.  Although  the  above  formula  is  approximate,  it  should  be 
sufficiently  accurate  for  purposes  of  the  present  study  since  changes  in 
drag  will  be  small  due  to  the  small  waterplane  area,  and  due  to  the  small 
changes  of  immersion  encountered. 
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A  vertical  force  model  was  developed  based  on  data  provided  by 
DTNSP.DC.  Due  to  the  wide  scatter  in  the  data,  the  only  variable  having  a 
discernable  effect  on  z  appears  to  be  the  speed.  An  analysis  of  the  data 
provided  leads  to  the  following  model: 

Z  =  (1.53  +  0.125U)pL2U2/2000. 

The  magnitude  of  the  above  force  results  in  only  minor  changes  in  heave. 

For  example,  at  a  speed  of  25  knots,  the  vertical  force  will  result  in  a 
change  in  heave  of  approximately  0.25  meters. 

A  hull  pitch  moment  is  estimated  from  the  drag  force  times  an 
appropriate  moment  arm.  The  moment  arm  is  chosen  such  that  the  resulting 
pitch  moment  is  zero  during  straight  ahead  flight  at  zero  degree  trim. 

RUDDER  FORCE  AND  MOMENT  MODEL 

Experimental  data  for  the  rudder  X,  Y,  K  and  N  forces  were  provided 
by  DTNSRDC  in  a  form  similar  to  the  hull  force  data.  Only  two  parameters 
were  varied,  the  rudder  angle  of  attack  and  the  speed.  Since  the  hull  data 
contained  the  effect  of  zero  rudder  angle,  the  rudder  force  described  in  this 
section  represents  the  increment  in  force  due  to  non-zero  angle  of  attack 
on  the  rudder.  From  geometric  considerations,  the  angle  of  attack  on  the 
rudder  used  in  the  rudder  force  model  equals  the  rudder  deflection  angle. 

If,  instead,  the  rudder  forces  had  not  been  included  in  the  hull,  the 
appropriate  rudder  angle  would  be  the  actual  inflow  angle  on  the  rudder, 
dependent  not  only  on  the  deflection  angle  but  also  on  the  ship  state  variables, 
such  as  drift  angle  and  yaw  rate. 
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For  various  rudder  sizes,  a  combination  of  theory  and  empirical 
considerations  yields  the  equation 


F  = 


Fo* 


Area ' 

X  Area^ 
o 


where  Fq  is  the  force  or  moment  for  the  design  rudder,  and  AR  =  aspect 
ratio.  From  pure  theory,  at  very  high  aspect  ratios  the  rudder  force  is 
independent  of  small  changes  in  the  aspect  ratio  and  is  proportional  only 
to  the  area  change.  At  very  low  aspect  ratios,  the  rudder  force  is  propor¬ 
tional  to  both  area  change  and  the  percent  aspect  ratio  change.  Since  the 
curve  is  smooth  between  these  two  points,  for  a  rudder  between  these  two 

/arAk 

extremes  the  rudder  force  is  proportional  to  l  AR~ i  where  K  lies  between 
zero  and  one.  Based  on  an  examination  of  experimental  data  and  theoretical 
considerations,  a  value  of  0.5  is  assigned  to  K.  Through  the  use  of  this 
formula,  changes  in  rudder  forces  due  to  pitch  and  immersion  changes  are  in¬ 
cluded  in  the  rudder  force  model.  The  rudder  Z  force  is  assumed  to  be  zero, 
and  the  pitch  moment  is  calculated  as  the  X  force  times  a  moment  arm  equal  to 
one  half  the  usbmerged  rudder  draft  plus  the  distance  from  the  coordinate 
system  origin  to  the  waterline.  The  above  formula  is  approximate.  However, 
only  small  changes  in  aspect  ratio  and  areas  are  considered  in  this  study. 
Under  these  conditions,  the  above  formula  should  be  sufficiently  accurate. 
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Data  was  also  supplied  for  the  fixed  forward  rudder  (Case  E  and  F). 

By  subtracting  the  Case  A  (no  fixed  forward  rudder)  forces  from  the  Case  E 
forces,  one  obtains  a  data  base  for  the  fixed  forward  rudder  as  a  function  of 
the  three  parameters:  speed,  drift  angle  and  yaw  rate.  The  drag  data  for  the 
fixed  rudder  are  quite  scattered.  Because  of  this,  the  theoretical  approxi¬ 
mations  relating  drag  to  lift  is  used:  drag  =  lift/7.  Since  the  additional 
drag  of  the  rudder  is  but  a  small  fraction  of  the  total  drag,  this  approxi¬ 
mation  is  adequate. 

A  three  dimensional  interpolation  procedure  is  used  to  obtain  the 
appropriate  non-dimensional  fixed  rudder  forces.  It  should  be  noted  that 
only  one  fixed  rudder  will  be  used  at  a  time,  the  starboard  one  for  starboard 
turns  and  vice  versa.  During  the  transition  time  of  zero  extension  to  full 
extension,  the  forces  are  assumed  to  be  in  direct  proportion  to  the  ratio  of 
extended  length/full  length. 

FIN  FORCE  AND  MOMENT  MODEL 

The  SWATH  configurations  considered  in  this  study  have  four  hori¬ 
zontal  fins  located  on  the  inside  of  the  hulls,  two  forward  and  two  aft. 

The  angle  of  attack  of  these  fins  can  be  individually  changed  to  control  pitch, 
induce  roll,  and  control  depth.  The  mathematical  model  used  to  characterize 
the  fin  lift  is  given  by, 

2 

Lift  =  (l/2)pAU/  ct 


where 


A  =  fin  area  or  chord  X  span 
Uf  =  local  inflow  velocity 
3C, 

— -=  lift  coefficient  slope, 
act 
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a  =  fin  angle  of  attack. 


The  value  of  the  lift  coefficient  slope  is  obtained  from  data  supplied  by 
DTNSROC.  The  values  for  the  slope  remain  constant  for  fin  angles  of  attack 
up  to  15°,  after  which  stall  occurs  resulting  in  a  rapid  decrease  in  lift. 

For  fin  angles  of  attack  less  than  15°,  there  is  good  correlation  between  the 
measured  value  of  1.19  and  the  theoretical  value  of  1.3.  Since  the  slope  of  the 
lift  coefficient  agrees  quite  well  with  the  theoretical  prediction,  it  is 
assumed  that  the  actual  lift/drag  ratio  can  be  adequately  approximated  by  the 
theoretical  prediction  for  similar  fins.  The  increment  in  drag  due  to 
non-zero  fin  angle  of  attack  is  obtained  directly  from  the  lift/drag  ratio 
and  the  lift.  The  fin  drag  at  zero  angle  of  attack  is  contained  in  the  hull 
force  data  base. 

The  fin  angle  of  attack  is  dependent  on  the  fin  deflection  angle 
and  the  ship  velocities.  The  local  velocities  are  given  by 

=  U  -  r.Yf  +  q-Zf, 

=  W  -  p-Zf  +  r-Xf, 

where  Xf,  Y^,  and  are  the  coordinate  location  of  the  fin.  It  is  assumed 
that  the  fluid  flow  parallel  to  the  fin  span  has  a  negligible  effect  on  the 
fin  lift  and  drag.  Since  the  fin  lift  and  drag  are  perpendicular  and  parallel 
to  the  stream  flow,  these  values  are  transformed  to  the  body  coordinate  system. 
The  force  along  the  body  y  axis  is  assumed  to  be  zero.  The  roll,  pitch,  and 
yaw  moments  are  obtained  from  the  transformed  forces  by  multiplying  by 
appropriate  moment  arms.  The  center  of  pressure  is  assumed  to  be  at  the  one 
quarter  chord  and  one  half  span. 

Note  that  the  fin  model  is  valid  only  up  to  angles  of  attack  of 
approximately  15°.  For  larger  angles  of  attack,  the  model  overestimates  the 


lift  and  under  estimates  the  drag  since  stall  effects  are  neglected.  During 
normal  maneuvers,  fin  angles  of  attack  larger  than  15°  are  rarely  required. 
The  results  of  work  being  conducted  by  Dr.  C.  M.  Lee*,  DTNSRDC,  will  be 
incorporated  at  a  future  date  prior  to  carrying  out  simulation  studies  using 
the  fins.  This  effort  presents  refined  estimates  appropriate  for  the  fins. 


C.  M.  Lee  and  Mr.  Martin,  "Determination  of  the  Size  of  Stabilizing  Fins  for 
Small  Waterplane  Area  Twin-Hull  Ships",  DTNSRDC  Report  4495,  Nov  1974. 
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PROPELLER  THRUST  MODEL 

The  SWATH  configurations  considered  in  this  study  are  driven  by 
two  submerged  variable  pitch  propellers.  The  mathematical  model  used  to 
represent  the  thrust  is  based  on  experimental  data  provided  by  OTNSRDC. 

These  data  were  obtained  for  a  single  SWATH  strut/hull /propeller  configuration. 
The  value  for  and  Kq  were  measured  for  advance  ratios  values  (J)  ranging 
from  -2.5  to  5.0  and  for  pitch/diameter  ratios  ranging  from  -4  to  4,  Data 
were  extracted  from  the  curves  and  are  used  in  a  two  dimen¬ 
sional  interpolation  procedure  to  obtain  values.  The  corresponding  thrust 
is  given  by  the  formula, 

2  4 

T  =  K_prTD  , 

where 


p  *  density  of  water, 
n  =  rotational  speed  of  the  propeller  in 
revolutions  per  second. 

0  =  propeller  diameter. 

The  thrust  is  assumed  to  act  along  the  thrust  line  used  in  the  experimental 
testing,  which  was  parallel  to  the  body  x  axis.  In  the  simulation  model, 
each  propeller  is  treated  separately  so  that  differential  thrust  can  be  applied. 
The  advance  coefficient  is  given  by 

J  -  Va/(nD), 

where  V.  is  the  speed  of  advance  of  the  propeller.  The  value  of  V,  is  less 

a  a 
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than  the  ship  velocity  due  to  hull  interference  effects.  These  hull  inter¬ 
ference  effects  are  taken  into  account  via  a  wake  fraction.  An  average 
value  of  0.82  is  used  in  the  mathematical  model,  i.  e., 


V  =  0.82  U. 

a 

The  propeller  forces  along  the  body  y  and  z  axes  and  the  roll  moment 
are  assumed  negligible  and  not  modeled.  The  pitch  and  yaw  moments  are  cal¬ 
culated  from  the  propeller  thrust  times  the  appropriate  moment  arms. 

The  propeller  model  neglects  changes  in  thrust  due  to  non-zero 
inflow  angle.  This  effect  is  small  during  high  speed  turns  but  may  become 
significant  for  low  speed  turns  that  result  in  large  drift  angles.  The 
large  changes  in  inflow  angle  can  lead  to  non-zero  values  for  the  propeller 
forces  along  the  y  and  z  body  axes.  The  effect  would  be  largest  for  low 
speed  turns  (less  than  5  knots)  carried  out  by  the  application  of  a  differ¬ 
ential  thrust.  In  these  cases  the  ship  transverse  velocity  can  be  a  signi¬ 
ficant  fraction  of  the  forward  speed.  Including  the  effects  of  changing 
inflow  angles  is  beyond  the  scope  of  the  present  effort.  It  is  expected 
that  these  effects  will  have  a  small  effect  on  the  predictions  presented  in 
this  report. 

AERODYNAMIC  FORCE  MODEL 

At  this  time,  the  forces  and  moments  due  to  aerodynamic  effects 
have  been  ignored  and  set  equal  to  zero.  A  rough  estimate  of  the  magnitude 
of  the  aerodynamic  forces  and  moments  indicates  that  a  25  knot  wind  at  90° 
would  result  in  a  roll  angle  change  of  less  than  0.25°.  Based  on  this  estimate, 
neglecting  the  aerodynamic  forces  and  moments  is  a  reasonable  approximation. 


22 


BUOYANCY  FORCE  AND  MOMENT  MODEL 


The  buoyancy  forces  and  moments  are  determined  from  coefficients 
Zg,  Zz,  Mg,  and  Mz  provided  by  DTNSRDC.  The  mathematical  expressions  for 
the  forces  and  moments  are  given  by, 

ZB  =  Ze0  +  ZzZ 
MB  =  M00  +  MzZ 

KD  =  (Beam/2 )2tan$  Z  . 

D  Z 

where  Zz  is  the  combined  value  for  both  hulls.  The  drag,  y  axis  force,  and 
yaw  moment  due  to  buoyancy  are  assumed  zero. 

RATE  DEPENDENT  FORCES  AND  MOMENTS 

Various  first  and  second  order  rate  dependent  coefficients  are 
contained  in  the  equations  of  motion.  These  are  of  the  form  X.  and  X.  .  where 

l  l  J 

X  represents  a  force  or  moment  and  i  and  j  represent  u,v,w,p,q,r.  The  entire 
term  is  the  local  derivative  of  X  with  respect  to  i(j)  so  the  associated  force 
or  moment  is  equal  to  X -  -  -Tj ,  where  T  and  j  represent  the  velocities  of  the 

1  J 

ship.  At  the  present  time,  experimental  values  have  not  been  provided  so 
all  terms  have  been  set  to  zero.  Research,  that  should  lead  to  theoretical 
estimates  of  these  coefficients  is  currently  underway  at  DTNSRDC  by  C.  M. 

Lee,  et  al . 
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SIMULATION  PREDICTIONS 


The  digital  simulation  system  developed  during  this  study  has  the 
capability  of  simulating  various  maneuvers  via  simulated  pilot  actions  or 
commands.  The  possible  pilot  actions  consist  of  rudder  deflections,  propeller 
RPM  changes,  propeller  pitch  changes,  fin  angle  changes,  and  changes  in  the 
position  of  the  forward  turning  foil.  All  controlled  elements  can  be  inde¬ 
pendently  moved;  i.  e. ,  the  starboard  RPM  can  be  set  to  0.0  and  the  port  RPM 
to  a  non-zero  value  simultaneously.  The  occurrence  of  a  pilot  action  is 
governed  by  the  value  of  time  and/or  state  variables.  The  pilot  actions  were 
used  to  carry  out  turning  maneuvers  at  high  and  low  speeds,  stopping  at  low 
speeds,  and  zigzag  type  of  maneuvers.  The  results  of  the  turning  and  stopping 
simulations  are  reported  in  this  section. 

Several  different  SWATH  ships  having  various  masses  and  slightly 
varying  length/beam  ratios  were  simulated.  Minor  modifications  in  the  mathe¬ 
matical  models  were  made  to  compensate  for  the  dissimilarities  between  the 
experimental  test  model  and  the  actual  ship  being  considered.  These  modifi¬ 
cations  were  discussed  in  detail  with  DTNSRDC  personnel  prior  to  generating 
final  results.  Results  for  one  set  of  ship  parameters,  as  listed  in  Table  5, 
are  presented  in  this  report. 
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TABLE  5  --  Design  Characteristics  of  Full  Scale  SWATH 
Used  for  Sample  Predictions 


Linear  ratio  full  scale/model 
Length 

Strut  Length 
Strut  Thickness 
Transverse  Strut  Separation 
Displacement,  design  draft 
Displacement,  deep  draft 
Design  draft 
Deep  draft 

Wetted  surface,  design  draft 
Wetted  surface,  deep  draft 
Longitudinal  distance,  nose  to  origin 
Vertical  distance,  baseline  to  origin 

Forward  fins,  each: 

Chord 

Span 

Aspect  ratio 

Longitudinal  distance,  origin  to  center  of  pressure 
Transverse  distance,  centerline  to  center  of  pressure 
Vertical  distance,  origin  to  center  of  pressure 
9C|  ,  per  radian 
3a 

Aft  fins,  each: 

Chord 

Span 

Aspect  ratio 

Longitudinal  distance,  origin  to  center  of  pressure 
Transverse  distance,  centerline  to  center  of  pressure 
Vertical  distance,  origin  to  center  of  pressure 
3Cl,  per  radian 
3a 


General 

Parameters 

22.5 

73.15m 

52.5m 

2.21m 

22.86m 

2.9x1 O^kg 

3. 01 5x1 O^kg 

8.13m 

9.27m 

2847m2 

3088m2 

35.14m 

10.36m 


2.61m 

2.903m 

1.11 

18.42m 

7.76m 

8.06m 

1.9 


4.  5m 

5.0m 

1.11 

26.71m 

7.76m 

8.06m 

1.19 
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TABLE  5  --  Design  Characteristics  of  Full  Scale  SWATH 
Used  for  Sample  Predictions  (cont.) 

General 

Parameters 


Propellers: 

Diameter  4.36m 

Nominal  Pitch/Diameter  ratio  1.0775 

Pitch/Diameter  range  -0.6  to  1.2 

Longitudinal  distance,  propeller  to  origin  39.72m 

Transverse  distance,  propeller  to  centerline  11.43m 

Vertical  distance,  propeller  to  origin  8.06m 

Strut  Rudder,  each: 

Chord  5.72m 

Span  7.32m 

Thickness  1.71m 

Span,  to  design  waterline  3.0m 

Area,  to  design  waterline  21.96m 

Longitudinal  distance,  origin  to  center  of  pressure  19.92m 

Transverse  distance,  centerline  to  center  of  pressure  11.43m 

Vertical  distance,  origin  to  center  of  pressure  --  3.73m 

design  waterline 

Spade  Rudder,  each: 

Chord  5.42m 

Span,  leading  edge  10.82m 

Span  to  design  waterline,  leading  edge  3.85m 

Thickness  0.813m 

Area,  to  design  waterline  21.9m 

longitudinal  distance,  origin  to  center  of  pressure  30.4m 

Transverse  distance,  centerline  to  center  of  pressure  11.43m 

Vertical  distance,  origin  to  center  of  pressure  --  4.1m 

design  waterline 

Fixed  forward  turning  foil: 

Chord  3.49m 

Span  6.98m 

Aspect  ratio  2.0 

Area  24.34m2 

Thickness  0.15 

Longitudinal  distance,  origin  to  center  of  pressure  18.42m 

Transverse  distance,  centerline  to  center  of  pressure  11.43m 

Vertical  distance,  origin  to  center  of  pressure  13.85 

Fixed  angle  of  attack  5° 


HIGH  SPEED  TURNS 


High  speed  turn  maneuvers  are  carried  out  by  deflecting  the  aft  rud¬ 
ders  from  zero  degrees  to  the  desired  rudder  angle  at  a  rate  of  5  deg/sec. 

The  simulation  is  run  for  a  sufficiently  long  time  such  that  the  heading 
angle  changes  by  180°. 

The  desired  turn  diameter  at  high  speeds  was  specified  to  be  less 
than  seven  ship  lengths.  A  comparison  of  turn  diameters  versus  rudder  angle 
for  strut  and  spade  rudder  configurations  indicates  that  the  specification 
can  be  met  by  the  spade  rudder  configurations,  but  not  by  the  strut  rudder  con 
figurations.  A  comparison  of  turning  performance  at  20  kts  for  both  configura 
tions  is  presented  in  Figure  5.  The  turn  diameter  specification  is  met  at 
speeds  of  25  kts  with  the  spade  rudder  configuration,  as  illustrated  in 
Figure  6 . 

If  turn  maneuvers  are  performed  at  a  draft  deeper  (1  meter)  than  the 
design  draft,  a  significant  reduction  in  turn  diameter  is  achieved  (see 
Figure  6).  However,  the  speed  reduction  during  the  turn  is  greater  for  the 
deeper  draft.  At  25  kts,  there  is  a  19  percent  reduction  in  speed  at  the 
deeper  draft  compared  with  a  12  percent  reduction  at  the  design  draft.  The 
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DIAMETERS  DIAMETERS  FOR  SPADE  RUDDERS 


percent  reduction  in  speed  is  greater  at  the  lower  speeds,  27  percent  for 
the  deeper  draft  and  16  percent  for  the  design  draft.  A  deeper  draft  can  be 
obtained  by  using  the  four  fins  to  induce  a  net  downward  force. 

An  analysis  of  Figures  5  and  6  indicates  that  there  is  a  small  in¬ 
crease  in  turn  diameter  as  the  rudder  angle  increases  beyond  approximately  35°. 
This  increase  in  turn  diameter  is  attributed  to  rudder  stall  effects.  Other 
factors  associated  with  the  fluid  flow  pattern  may  also  contribute  to  this 
effect.  i 


The  excursions  in  roll  angle  and  drift  angle  are  minor  in  all  cases 
simulated.  For  example,  the  steady  state  roll  angle  in  a  20  Kt  turn  is  less 
than  4°  and  the  drift  angle  is  less  than  2.5°.  Mote  that  the  ship  rolls  out 
during  a  turn. 


The  forward  turn  foils  are  not  employed  in  the  above  simulation 
studies.  The  use  of  these  rudders  may  induce  a  sufficient  turning  moment  so 
that  the  strut  rudders  would  be  sufficient  to  turn  the  craft  within  seven 
boat  lengths.  A  suitable  mathematical  model  necessary  to  model  these  forward 
rudders  or  foils  is  currently  under  development.  It  is  expected  that  the 
relative  differences  between  the  curves  appearing  in  Figures  5  and  6  will  re¬ 
main  the  same  when  the  forward  turn  foils  are  employed.  The  net  effect  should 
be  to  lower  all  of  these  curves  by  approximately  the  same  amount. 

A  comparison  of  SWATH  turn  trajectories  at  20  kts  with  a  destroyer 
turn  trajectory  at  approximately  the  same  speed  is  illustrated  in  Figure  7. 

The  SWATH  transfer  is  approximately  25  percent  greater  than  the  transfer  for 
the  destroy er.  The  SWATH  turn  maneuver  employed  only  the  spade  rudders  at 
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CTORY  COMI’ARFI)  WITH  A  DISIROYIR  TURN  IRAJICTORY 


design  draft.  The  forward  turn  foils  or  differential  thrust  capabilities  are 
not  utilized. 

The  SWATH  turning  performance  at  low  speeds  is  quite  impressive. 

Low  speed  turn  results  are  discussed  in  the  next  subsection. 
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LOW  SPEED  TURNS 


Low  speed  turns  (5  and  10  kts)  are  accomplished  by  applying  a  dif¬ 
ferential  thrust  in  addition  to  using  the  aft  rudders.  The  objective  of  the 
simulation  effort  is  to  determine  the  proper  combination  of  differential  thrust 
and  rudder  angle  such  that  the  advance  is  a  minimum  for  a  90°  turn.  Also,  the 
horse  power  requirements  are  to  remain  less  than  a  specified  value.  The  de¬ 
sired  performance  can  be  accomplished  by  changing  the  propeller  pitch/diameter 
ration  on  one  side  of  the  ship  to  -0.6  and  the  opposite  pitch/diameter  ratio 
to  1.0.  The  resulting  predicted  advance  is  approximately  110  meters  for  a 
5  kt  turn.  At  the  lower  speeds,  the  rudders  are  not  very  effective  in  turning 
the  ship.  In  fact,  the  predicted  turn  diameter  achieved  by  applying  only  a 
differential  trust  is  about  15  meters  larger  at  5  kts. 

During  the  low  speed  turns,  both  the  drift  angle  and  non-dimensional 
yaw  rate  exceed  the  range  of  the  available  data  base.  The  drift  angle  may  in¬ 
crease  to  10°  and  the  non-dimensional  yaw  rate  to  values  up  to  1.5.  The 
reason  for  such  large  values  is  that  the  ship  is  turning  at  a  high  rate  at  a 
very  low  speed  and  is,  in  effect,  pivoting  about  the  propeller  location.  Ex¬ 
trapolation  of  the  data  should  provide  reasonable  results  due  to  the  extreme 
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linearity  of  the  data  at  low  speeds.  Further  experimental  tests  would  be 
necessary  to  fully  determine  the  validity  of  the  extrapolations. 

The  ship  advance  during  low  speed  turns  is  only  slightly  larger  than 
the  ship  length  (75m).  Hence,  in  obstacle  avoidance  maneuvers,  it  is  of  in¬ 
terest  to  examine  the  predicted  path  swept  out  by  SWATH,  in  addition  to  the  CG 
trajectory.  The  actual  path  swept  out  by  SWATH  during  a  low  speed  turn  is 
illustrated  in  Figure  8.  These  results  indicate  that  SWATH  is  capable  of 
avoiding  a  point  obstacle  approximately  112  meters  in  front  of  the  craft,  by 
a  turning  maneuver.  A  better  maneuver  for  avoiding  point  obstacles  at  low 
speeds  involves  stopping  the  ship  Stopping  maneuvers  are  discussed  in  the 
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STOPPING 


Stopping  maneuvers  were  performed  at  5  and  10  kts  by  changing  the 
propeller  pitch/diameter  ratio  to  a  value  of  -0.6*  over  a  period  of  5  seconds 
Without  deflecting  the  rudders  to  increase  drag,  the  predicted  stopping 
distances  are  27m  and  80m  for  speeds  of  5  and  10  kts  respectively.  However, 
reducing  the  pitch/diameter  ratio  to  a  value  of  -0.6  over  a  5  second  period 
results  in  higher  horsepower  requirements  then  desired.  In  order  to  maintain 
horsepower  requirements,  the  minimum  pitch  value  is  set  to  be  -0.15.  The 
resulting  stopping  distances  for  5  and  10  kts  are  34m  and  97m,  respectively. 

In  order  to  improve  stopping  performance,  the  rudders  can  be 
deflected  to  their  maximum  values  thereby  increasing  the  drag.  However,  the 
net  change  in  stopping  distance  is  negligible.  Such  a  maneuver  may  be  more 
effective  for  high  speed  stopping. 
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FIN  EFFECTS  ON  TURNS 


The  fins  can  be  employed  to  change  the  roll  angle  during  turn.  It 
was  found  that  by  forcing  the  ship  to  roll  out  more  during  a  turn,  the  turning 
rate  increased  slightly,  and  vice  versa.  Because  of  stall  effects,  the  fin 
angles  were  limited  to  values  less  than  15°  during  this  investigation.  At 
these  angles,  the  fins  have  very  little  effect  on  the  turn  diameters.  Based 
on  the  comparisons  of  predicted  turn  diameters  for  different  drafts,  it  may  be 
more  beneficial  to  use  the  fins  to  produce  a  net  downward  force  and  thus  in¬ 
crease  the  draft. 


DIGITAL  SIMULATION  DESCRIPTION 


OVERVIEW 

The  digital  simulation  system  used  to  implement  the  mathematical 
model  for  numerical  solution  is  a  FORTRAN  computer  program  designed  to 
operate  on  the  CDC  6700  digital  computer  located  at  ITNSRDC,  Carderock.  The 
final  code  is  based  on  a  preliminary  design  developed  using  a  "top  down" 
approach.  Many  features  of  existing  digital  simulation  programs  for  AALC  are 
incorporated  in  the  present  system,  such  as,  the  pilot  action  strategy  (Sub¬ 
routine  CONTRL),  hydrodynamic  table  lookup  procedure  (Subroutines  HYDINT, 
HYDPAR) ,  and  the  Runge-Kutta-Merson  integration  algorithm  (Subroutine  KUTMER)t 
Several  of  the  modules  employed  have  been  previously  developed  and  tested  by 
ORI  (TAB2V,  TAB3V,  NGHBR1) .  The  highly  modular  structure  of  the  digital  simu¬ 
lation  system  with  well  defined  interfaces  allowed  simultaneous  and  on  sched¬ 
ule  development  and  unit  testing  of  all  force  and  moment  subroutines  by 
several  different  analysts.  Also,  the  modular  design  is  extremely  flexible 
allowing  various  SWATH  configurations  to  be  easily  modeled. 

A  block  diagram  for  the  complete  simulation  system  is  illustrated  in 
Figure  9.  The  order  in  which  programs  are  called  is  in  a  downward  fashion 
from  left  to  right,  except  the  block  data  module  which  is  initialized  at 
compile  time.  The  main  program  serves  as  a  driver  for  the  system  as  illus¬ 
trated  by  the  Flowchart  appearing  in  Figure  10.  Brief  descriptions  of  major 
subroutines  and  functions  are  presented  in  Table  6.  More  detailed  descrip¬ 
tions  of  the  subroutines  are  presented  later  in  this  report.  A  listing  is 
presented  in  Appendix  B. 

*R.F.  Messalle, "First  Generation  Dynamic  Simulation  Model  for  an  Air  Cushion 
Vehicle",  ORI,  Report  #814,  September  28,  1973,  Contract  No.  N00167-73-C-0213 


MAP 


fiyure  9  -  SWATH  Digital  Simulation  System  Block  Diagram 


ENTER 


1 1 


END 


Figure  10  -  Flowchart  for  Main  Program 
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TABLE  6 

SWATH 

INPUT 

OUTPUT 

CONTRL 

TURNON 

DAUX 

BLOCKD 

ADDFOR 

HYRITE 

DRAAGG 

HULFOR 

BOYFOR 

FINFOR 


-  SUMMARY  OF  SUBROUTINE  DESCRIPTIONS 

Main  program.  Serves  as  a  driver  for  the  entire  system. 

Reads  control  data,  NAMELIST  input  and  hydrodynamic  data, 
converts  variables  to  internal  MKS  units  and  initializes 
state  variables. 

Loads  output  data  into  page  arrays  and  writes  results  on 
line  printer. 

Control  system.  Changes  control  surfaces  according  to  pilot 
action  data. 

Turns  on  all  intermediate  printout  options  used  for  de¬ 
bugging. 

Computes  mass  matrix,  calls  all  force  subroutines  and  com¬ 
putes  total  force  vector.  Calls  LINQ  to  solve  simultaneous 
equations  for  state  acceleration  vector. 

Block  data.  Initializes  all  pertinent  parameters.  Default 
values  can  be  overidden  by  processing  in  subroutine  INPUT. 
Computer  forces  and  moments  due  to  rate  dependent  terms  in 
operations  of  motion  (Nr,  Yy,  etc.). 

Prints  out  a  table  of  the  hull  hydrodynamic  Y,  K  and  N  forces 
and  moments. 

Computes  the  drag  (X  force)  on  the  hull. 

Computes  the  hull  hydrodynamic  forces  and  moments  by  calling 
appropriate  subroutines. 

Computes  forces  and  moments  due  to  buoyancy. 

Computes  all  forces  and  moments  due  to  the  fins. 
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TABLE  6 
TOTVEL 

RUDFOR 

AIRFOR 

PROFOR 

LINQ 

EULER 

TAB2V 

NGHBR1 

TAB3V 

HYDINT 

HYDPAR 

KUTMER 


I 

-  SUMMARY  OF  SUBROUTINE  DESCRIPTIONS  (continued) 

Calculates  total  inflow  U,  V,  W  velocities  given  the  position 
vector  and  the  ship  motions. 

Calculates  forces  and  moments  due  to  aft  movable  rudders  and 
fixed  forward  rudder. 

Calculates  aerodynamic  forces  and  moments. 

Calculates  propeller  forces  and  moments. 

Solves  for  X*  where  tx  =  t.  In  this  case,  t  =  (u,  v,  w,  p, 
q,  r )• 

Performs  Euler  integration  on  state  vector  and  updates  time. 
Subroutine  to  perform  2-dimensional  interpolation  on  data. 
Finds  nearest  neighboring  points  in  grid  arrays  used  by  inter¬ 
polation  subroutines  TAB2V  and  TAB3V. 

Subroutine  to  perform  3-dimensional  interpolation  on  data. 
4-way  interpolation  subroutine  to  calculate  hydrodynamic  hull 
forces  based  on  speed,  roll  angle,  drift  angle  and  yaw  rate. 
Computes  interpolation  coefficients  for  HYDINT. 

Subroutine  to  perform  Runge-Kutta-Merson  integration. 
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Conmuni cation  between  various  subroutines  is  handled  via  named 
common  blocks  wherever  possible.  The  FORTRAN  symbols  appearing  in  the  common 
blocks  are  defined  in  a  global  sense.  For  example,  the  symbol  U  represents 
the  ship  velocity  along  the  body  x-axis  in  all  subroutines.  Definitions  of 
all  variables  appearing  in  named  common  blocks  are  presented  in  Table  7.  The 
table  contains  brief  descriptions,  the  type  of  variable  (Real,  integer,  etc.), 
and  the  units.  Generally,  descriptive  notation  is  employed  in  defining  the 
FORTRAN  symbols.  Some  exceptions  are  (XC,  YC,  ZC)  which  represent  (XG,  YG, 
ZG),  the  CG  location. 

Core  requirements  are  approximately  41,000  words  base  10.  Execution 
time  is  generally  5  times  or  more  faster  than  real  time. 

A  definition  of  required  job  control  commands,  a  definition  of  all 
input  data,  and  a  user's  guide  are  presented  in  the  next  section.  Several 
test  run  setups  are  also  presented. 
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TABLE  7  DEFINITION  OF  FORTRAN  SYMBOLS  APPEARING  IN  NAMED  COMMON  BLOCKS 


FORTRAN 

Symbol 

Descri otion 

mm 

Uni  ts 

AIRFM  (1-3) 

COMMON  /AIRVAR  / 

Air  resistance  X,  Y,  Z  forces  along  body  axes 

Real 

nts 

AIRFM  (4-6) 

Air  K,  M,  N  moments  about  body  axes 

Real 

n;-m 

ADMTRX(6,6,6) 

COMMON  /ADDVAR  / 

Matrix  of  rate  dependent  2nd  order  terms  for 

Real 

ADDFM  (1-3) 

hull  forces 

Rate  dependent  X,  Y,  Z  forces 

Real 

nts 

ADDFM  (4-6) 

Rate  dependent  K,  M,  N  moments 

Real 

nt-m 

RTMTRX(6,6) 

Matrix  of  1st  order  rate  dependent  terms  for 

Real 

- 

MSIIBZ 

hull  forces 

COMMON  /BOYVAR  / 

Restoring  derivative  Mz  due  to  buoyancy 

Real 

nt-m/m 

MSUBTH 

Restoring  derivative  MQ  due  to  buoyancy 

Real 

nt-m/rad 

ZSUBZ 

Restoring  derivative  Zz  due  to  buoyancy 

Real 

nts/m 

ZSUBTH 

Restoring  derivative  Zg  due  to  buoyancy 

Real 

nts/rad 

BUOYFM  (1-3) 

Buoyancy  X,  Y,  Z  forces  along  body  axes 

Real 

nts 

BUOYFM  (4-6) 

Buoyancy  K,  M,  N  moments  about  body  axes 

Real 

nt-m 

XC,YC,ZC 

COMMON  /CGMOVE  / 

X,  Y,  Z  coordinates  of  center  of  gravity 

Real 

m 

SC FACT 

COMMON  /CRAFT  / 

Ship-Model  scale  factor 

Real 

Ship  length 

Real 

m 

MASS 

Ship  mass 

Real 

kg  2 

IXX, IYY.IZZ 

Ship  moments  of  inertia 

Real 

kg-mg 

IXY ,  IYZ.IYZ 

Ship  products  of  inertia 

Real 

kg-ni 

RADIAN 

COMMON  /CONST  / 

Number  of  degrees  per  radian,  57.29 

Real 

deg  3 
kg/m^ 

RHOA 

Density  of  air 

Real 

RHOA2 

4  of  RHOA 

Real 

m  ifiuf 

K  jffl 1 

RHOW 

Density  of  water 

Real 

RHOW2 

4  of  RHOW 

Real 

FTTOM 

Conversion  feet  to  meters 

Real 

m/ft 

TABLE  7  (Continued) 


Description 


COMMON  /OOTCQE  / 

Dimensional  coefficient 

Dimensional  coefficient  Y^ 

Dimensional  coefficient  Y- 

P 

Dimensional  coefficient  Y- 


Dimensional 

Dimensional 

Dimensional 


coefficient  Z- 

to 

coefficient 
coefficient  K- 


Dimensional  coefficient  K- 


Di mensional 
Dimensional 


coefficient  K* 
coefficient  M- 


Dimensional  coefficient  M- 


Dimensional 

Dimensional 


coefficient  N; 

v 

coefficient  N; 


Dimensional  coefficient  N- 

r 

COMMON  /FINVAR  / 

X,  Y,  Z  coordinates  of  fin  center  of 

pressure  1  =  bow  starboard,  2  =  bow  port, 

3  =  aft  starboard,  4  =  aft  port 

Areas  of  fins 

Lift  coefficient  of  fins 

Lift/drag  ratio  of  fins 

Fin  angles 

Fin  X,  Y,  Z  forces  in  body 

Fin  lift  and  drag 

Fin  X,  Y,  Z  forces  in  body  system 
Fin  K,  M,  N  moments 
Fin  angles  at  time=0 


TABLE  7  (continued) 


FORTRAN 

Symbol 

Description 

Type 

Units 

^  HULLFM  (1-3) 

COMMON  /HULVAR  / 

Hull  X,  Y,  Z  forces  along  body  axes 

Real 

nts 

HULLFM  (4-6) 

Hull  K,  M,  N  moments  about  body  axes 

Real 

nt-m 

HULPIK 

Difference  in  draft  from  design  condition. 

Real 

m 

LG  (4) 

deeper  draft=positi ve 

COMMON  /HYDRO  / 

Number  of  grid  points  for  HYDINT 

Integer 

LZ 

Number  of  zero  speed  curves 

Integer 

HY  (4,3213) 

Data  points  for  HYDINT  array 

Real 

nt  or 

XG  (4,  20) 

Grid  points  for  HYDINT  array 

Rea": 

nt-m 

m/s,  rad 

I POINT  (4) 

Pointer  for  HYDINT 

Integer 

- 

OT 

COMMON  /INTGRL  / 

Step  size  for  integration  orocedure 

Real 

s 

DTMIN 

Minimum  step  size  for  KUTMER 

Real 

s 

TMAX 

Maximum  time  for  run 

Real 

s 

AE  (15) 

Absolute  error  criterion  for  KUTMER 

Real 

- 

EPS  (15) 

integration 

Relative  error  criterion  for  KUTMER 

Real 

_ 

NSTATE 

integration 

Number  of  state  variables 

Integer 

- 

|  I DOF 

COMMON  /OPTION  / 

Degrees  of  freedom:  4,  5,  6 

Integer 

! IOPTN ( 15 ) 

IOPTN  (2)  3  0,  Do  not  print  pages  3-5 

Integer 

- 

t 

i 

PDIA 

IOPTN  (2)  =  1,  Print  pages  3-5;  rudder, 
fin,  and  propeller  forces 

COMMON  /PROPVR/ 

Prop  diameter 

Real 

m 

PROARM  (3) 

Prop  X,  Y,  Z  coordinates 

Real 

m 

PRPMS 

Prop  rpm  starboard 

Real 

rpm 

PRPMP 

Prop  rpm  port 

Real 

rpm 

PITCHS 

Pitch/Diameter  ratio  -  starboard 

Real 

- 

PITCHP 

Pitch/Diameter  ratio  -  port 

Real 

- 

PROPFM  (1-3) 

Prop  X,  Y,  Z  forces  along  body  axes 

Real 

nts 

PROPFM  (4-6) 

Prop  K,  M,  N  moments  about  body  axes 

Real 

nt-m 

PTCHSO 

P/D  ratio  -  starboard  Time  =  0 

Real 

- 

PTCHPO 

P/D  ratio  -  port  Time  =  0 

Real 
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TABLE  7  (Continued) 


4 


FORTRAN 

Symbol 

Description 

Type 

Units 

IR  (1-5) 

COMMON  /RDPRNT  / 

Read  units 

Integer 

IP  (1-10) 

Print  units,  6  =  main  output  printer 

Integer 

PNT  (15) 

Optional -debugging  output,  9  =  hydro  array 

Logical 

- 

RUDANG  (2) 

10  =  intermediate  force-moment  time  histories 

COMMON  /RUDVAR  / 

Aft  rudder  angles,  1  =  starboard,  2  =  port 

Real 

rad 

UUVEL  (10) 

Grid  velocity  array  for  rudder  forces 

Real 

m/s 

RALFA  (10) 

Grid  rudder  angle  array  for  rudder  forces 

Real 

deg 

AXRUD  (100) 

X  rudder  force  array 

Real 

nts 

AYRUD  (100) 

Y  rudder  force  array 

Real 

nts 

AKRUD  (100) 

K  rudder  moment  array 

Real 

nt-m 

ANRUD  (100) 

N  rudder  moment  array 

Real 

nt-m 

ARATIO 

Rudder  sizing  ratio,  to  standard  rudder 

Real 

- 

RUOANO  (2) 

Rudder  angles  at  time=0 

Real 

rad 

NUMUUV 

Number  of  data  grid  velocity  points 

Integer 

- 

NUMRAL 

Number  of  data  grid  rudder  angle  points 

Integer 

- 

BRUDAS 

Starboard  bow  rudder  position 

Real 

m 

BRUDAP 

Port  bow  rudder  position 

Real 

m 

RUDFM  (1-3) 

Rudder  X,  Y,  Z  forces  in  body  system 

Real 

nts 

RUDFM  (4,6) 

Rudder  K,  M,  N  moments  about  body  axes 

Real 

nt-m 

XRUD  (2) 

X  coordinates  of  aft  rudder  center  of 

Real 

m 

YRUD  (2) 

pressure 

Y  coordinates  of  aft  rudder  center  of 

Real 

m 

ZRUD  (2) 

pressure 

Z  coordinates  of  aft  rudder  center  of 

Real 

m 

RUDARM 

pressure 

Aft  rudder  pitch  moment  arm 

Real 

m 

BRDAPM 

Bow  rudder  pitch  moment  arm 

Real 

m 

BRSPAN 

Maximum  bow  rudder  extended  length 

Real 

m 

u.v.w 

COMMON  /STATE  / 

Ship  velocity  components  in  body  coordinate 

Real 

m/s 

P.Q.R 

system 

Sh i d  anaular  velocity  about  body  axes 

Real 

rad/s 

PHI 

Roll  angle 

Real 

rad 

THETA 

Pitch  angle 

Real 

rad 

PS  I 

Yaw  angle 

Real 

rad 

XE,YE,ZE 

Ship  position  in  fixed  space  coordinates 

Real 

m 

S 

Distance  travelled 

Real 

m 

SPEED 

Ship  speed 

Real 

rad 

DR I  FAN 

Drift  angle 

Real 

rad 
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TABLE  7  (continued) 


FORTRAN 


Symbol 


Description 


COMMON  /STATE  /  (Continued) 


COSPHI  ) 
COSTHE 
COSPSI  5 
SINPHI  i 
SINTHE 
SIMPS  I  ». 
TMTRX  (3,3) 

TIME 


Cosines  of  PHI,  THETA,  PS  I 


Sines  of  PHI,  THETA,  PSI 

Transformation  matrix  from  fixed  system  to 

body  system 

Time 


COMMON  /STATEO  / 


TO 
UO 
VO 
WO 
PO 

QO  / 
RO 

PHIO 
THETAO 
PS  10 


Initial  time 


Initial  values  of  u,  v  w,  P,  Q,  R,  0,  0,  w 


\ 


USER'S  GUIDE 


Examples  of  input  run  streams  for  Cases  A  -  D  are  illustrated  in 
Tables  8  -  11,  respectively.  Various  temporary  and  permanent  files  are 
required.  A  summary  of  the  files  follows: 

TAPE  8  ) 

TAPE  10  Temporary  files  to  store  intermediate  output 
TAPE  12  ) 

BINAR  -  Binary  file  for  the  system 
TAPE  1  -  Data  array  for  hull  hydrodynamic  Y  force 

TAPE  2  -  Data  array  for  hull  hydrodynamic  K  moment 

TAPE  3  -  Data  array  for  hull  hydrodynamic  N  moment 

HYDINT  -  Binary  file  continainq  subroutines  HYDINT 
and  HYDPAR 

The  asterisks  appearing  in  the  card  images  should  be  replaced  by 
account  numbers  and  ID  designator  appropriate  for  the  user.  Files  will  be 
copied  to  appropriate  ID  files  upon  request  through  Code  1572,  DTNSRDC. 

Four  types  of  data  are  read  from  cards:  a  descriptive  title, 
Namelist/SWATHD/,  pilot  action  data,  and  rate  dependent  coefficients.  De¬ 
scriptions  of  the  data  sets  follow. 

Title  Card 

The  first  data  card  is  a  title  card.  A  full  80  columns  may  be  used. 
Name!  i  s  t/ SWATH  D/ Da  ta 

Program  control,  except  for  pilot  actions,  are  governed  by  data 
located  in  Mamel i s t/ SWATH D/.  Namelist  input  is  used  to  minimize  data  set  up 
errors,  and  to  relate  directly  data  values  with  FORTRAN  symbols.  Many  of  the 
FORTRAN  symbols  are  initialized  in  the  BLOCK  DATA  program.  The  data  cards 
merely  override  the  values.  A  complete  description  of  all  symbols  appearing 
in  Namelist/SWATHD/  is  presented  in  Table  12.  Default  values  are  set  up  for 
ship  parameters  presented  in  Table  5.  A  few  points  should  be  noted  regarding 


TABLE  8 


-  SAMPLE  RUN  STREAM  FOR  CASE  A 
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DRAFT 


A  K  c  1 1  ■ '  r  '  *0.1. 

•.*’j  j  "  «c  ,  ^  A  ,5  .'i:  »or  «  -<*.  ~bt  + .  c  • 

- 1;. .  -» ;  r  +  ^  ~  . 

?»llr*0 -J  C  = -3  c.»  ’  .  ;  7L'  ♦  j '  «-2  7f  ^ 

-T.!lr*Ot, 

'*  •  0  a  Z  *  P  h  *  4  *  S  f  !I  +  7  t  .  4 . 4  1  £  +  ' .  c  ,  ~  .  c  r  i  ■*  n  ' ’  .  -  ~  %  r  n  l*  +  -  /- •  # 

-4.  t  r.~*5  &  , 

4»a'ir  +  0o%4«‘TCr_+,3''-;*4  .4  ?  r  +  ~  £, «  7  .3  r  •'  ♦  ?  .-.  , -  ■z  %  •*  t*:'  ,<  t 

-4.  nic+jc, 

—  O»40£+J  S.  . 

:.5  =  "+3A,5.f';-:'*  ,-4,'7a  +  '  av 

-5«f9“+Qfi,» 

3')R(J7:  v+n.0, 

4.02E+CF 
-4.  C2F.  +  0F,  * 

1  .Q  +  »  1  .  7*r-f  07,l."A-<"’7, 

-1  .  ”  ?r-*>0  7  , 

2»13F.+  .n7,2»',!*!r'»  )7»1  .-P~+'7»l  «^7‘.*i.|7«-l,r'Tr  +  "-',. 
-2.  1  '■=>!)  7  » 

2.4f>''*07.2.2^-.*l'7,2.C7r*}',,l.=  7r*C/,-l*37r  +  r7.. 

-? • 4nF +n  7 « 

-l.^Fp+77t 

4  *1  .  9  3F  *  :  7. 4  ♦-!  .  -  5~.  *  J  7 , 

'•My  =  12G  .  .c  Tru.p  '-i.  -•»  , 

r®P'*SO  =  l  *>.7  «  .  F3  1  -p?  rl  =0  .  ,  1 1- 

2  PILOT  i C T  T  ~  N  *3 

16  7  -7=.7 

i c  :/.°s  p  c .  r 

1300. r+oo 


_  ■  •  cor  •  o  ^  ^  a  /  »*  .  #  , 


-4.42r*Qf  ,  -4. 7-T  +  r 


-4 .  ■«  7  *3  4 .  -P.4  7 7«-  7e 


-r.'r:r*Qt,_S.C-.r4.C£ 


-l.Sfr*:7.-!.7t'*f7 


-,."s:‘,’7.-\:t'K' 


_  **  C77  .  ''I  <*■  <>  — 

J.  •  '  +  J  ‘  fc  ♦  '  / 


J  A T «-•  TrS  T 
t  S'4  a  7H7  VR'J^s-ft. 

7  p  T  LnT  fiCTTO1 

loo : .f+o  " 

1CQ',C0.G 

•r  OF  c’lL17 


TABLE  8  -  SAMPLE  RUN  STREAM  FOR  CASE  A  (CONTINUED) 
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j  =  c  *  '•  r. ,  c  i  r>  <;• «  c  « T  2"  a  .  r-  . 

Crft*c,r  '*’5*%  ♦**■*■* 

3~AU-"ST,T1?F1  1  »*°F. 

'>rCU",?T.--*  Fc  1  ',»r  r. 

-  -"Q  u  F  F T  »T  *3  "t  . 

iTTAC^-^iPLl  ,L^uvr>4T‘>  .*Y;t.rr 

47TAFFfTAPF?  .LA^vF/i  -f.  .-v-c,  r-._ 

4  tt  4CH,  ^APF^^Li  ■:-  y  Tr-  »'*Y  =  it  I  "  = 
\TTAO*- I  »  ,,!LJ  ?'•  *T~  "■  T  T  "  = 
ATTiC'-thY  FI'JT.--  T"  ■.?  v  *  {.  -  .  I  '  = 

Li4D ( HYC 1*7) 
u  r*;  as  . 

r  ar  v  *  T  ir’r  1  0  «  '  u  r«"  i  i  r  >  . 

F  - '  J  I  \!  ?  <  T  4  F  F 1 1  >  . 

C  '  -  Y  *  T  AFr  11  , 

° T  <i  i '  •  ^  <  r  a  p  r > . 

c;  py  <  tape®  ,cutpi.it>  . 

r '|  r\  -vc  c  p  r 

S ‘f a  T -J  TC3T  ~  Z'~  C 

*>’.</  i  TM2  '"•!>'!  =!♦  r^-r6.  F.5T'M 

AX"  1!^=  °*Z  . 


-  c  •  '1>=C-  PBTV(2l  =1,  ^ 


7  4!+f  A  ,-c  .3-  :  *3  a  .-A.  '  7  4  ,  -  4  .  7  7"  +  !  A  ,  a,  ;:•  7~  +  n  a  ,f  .  7.r+-;  4  ,  s 

}»74F-*-04* 


I.??7*®",?.!  -r*05.? 


'5.a"r-*-0='» 

-A.  3  7r*r  5  .-c  .  !0r-*-f.»5  7. .C',il*-OF  •?.  C  7 1  ♦  3  =  .  I.  o7E>"  r  «F 

-‘i.T'F+o* -6-  +  aF  .f  3r  =  ,r  +  r,c.^r*or  .7 

7.3 "r*CF  . 

-9.7'H>':  F  ,-7  ♦  -=.  "~£>  ?  =  .  . FFc*f'c  4  FF-*-?  "  .  =  .  7,1  F*?  c  ,7 

".74"*05  , 

-9.  44P+Q  =  ,-f, *C5  .-3."  •  77T*'-?.?*77‘>5-*,!.!5e£*CR,e 

9,44r  +  0c  « 

A  Y  0  !  J  2  =  o,o.~  . 

-1.  r.  4  F+0  c  »-l  .F7  - *0C  .-!.=  •  '♦  0=  ,-!  -+75 , 1.2  •'  •"♦P  c  *1.  *''■?*<'  c.l 

l.»4'>05  * 

-  6  •  F  !  E  ♦  u  F  •  -  F  »7"f.qc  . ;*  Lr+Cq»r.Alr+3  -  •‘:.c4r+f,c.'' 

(:  .5 1 r  *0F  » 

-1.  ?1E*CF  ,-1.3  7E*OX  ,-l.7  •  £+n-  z*-f-  .7At*2  =  .F.7i3f-»,,c.t.?lE*ft  *»1 

!.22r^0F  , 

-1.F7F4-0F.-1  l.r'c“*'  -  ,-l.*4r*?F.  1. 4i-^  '.  l.atr+o  e. 

l.F7r*00 . 

-1  .  '•  .-!  ."-F*C6,-1.7"F»  C4  .  -1  .1  17  *■  .  3.  lr  f :*;1  c  ,1.77F*f*A,? 

-l.Fs^:t,.1(<;:"40'  •- 1 . 7  .llr*:-F  .1.1’  .l.7  7F«-7  4,1 

ua4r*;F . 


TABLE  9  -  SAMPLE  RUN  STREAM  FOR  CASE  B 


DRAFT 


i  ,jn_  r.  »n  ,  1  . 

i  .0  4E>'v-. « '■•.?-£>  o^,  7.  y7z  ^7r*’ 

-i.  . « . 

+  .7.r  f «- ’r+"-. 

-  7.  o:  <■ , 

*  •  6  ‘i ” ♦  C ^  1  'r  *'f  *  *  T*  *Oft  Tf  +  0  1  . 

-  S  .  o  .  T 

7:  "♦Or  ♦c-.  T/1  r+~  «  1. 7i“*  ’ft  «  7.0  *■*"  «■  .-7.  -  -T*  :  r  ,-:.74r  +  :c,  _e.  7  i '+ C  -  . 

-8.  7'>t-  +  C  ?■  . 

7r  +  r,ft  .ft  .& -  -■»nft  .-ft . ft  -  •“«■  n r » '»•-'.  ?  ?r4. -ic , 

4  *rj  •  0  uE*  0  •'  *4  *-ft.  .  T  El  *  "  '  . 

*\°LtM  =  ^  »  C  •  ?  t 

7  *4 .2  7"  3. 11"  vf’ft  .-3.  1  1"«-  1  ft  *  3 ‘-ft.?7  -  *^ft. 

i. '.?r+')T ,?.n:r.>o 7,1  i  .c  ="  *o 7,-i ,”E>o'r •  .?r^*n7f ~r  +  r  -t, 

-l.r-?ft  4-0  7, 

*.T  !C*C7  » ■,.r'  7,T  .  :  ft  7*  "  .  -  7«-r>.  r~*l  7 ,  _  : .  7  r  +  7  7 ,  -  * .  =  -  -♦<•7, 

-7.  :  :r+j7. 

ft  .67"*C7 , 7.7;"*:  7.f  ,-4. 1  4F  +  ;  7.-f  .  -  xr  *n  7t  _7.  ->  .-r^f  7. 

-ft. ft  *F>0  7, 

<*5  ."  0r*  "7.  T.  ?'*■->  0  7  T.  ^ '*;*?  7,  7  r  +  .:.  7, 

**4 . 7ftE.*  C  7  •  7  .c4r*f>7  ft.  ;  -i!;+  7  7  7  7  »-«  .?  ft  7t 

T^Axsl^c  .  .->rr-P  -  =  i.C7.  n'-mir-j.r?, 

-  3  o  •  3  C  =  1 c  0  •  •  ocn-c3=]c3#,  it:''.’’.  ^yTsj  t  -ft*-. 

9  ft  I  LOT  iCTT,M3 

to  17.  ft  ft  ft  -ft  ft.  ft  ft,.  1 

1-  1  7 .  Q  U  a  -  ft  ft  .  ft  -  ,  ft  1 

-100 c  .r*o r 

- : :  0 "  7 : .  0 
3*‘-T^  TftST  r».ft 
“TV  t  ’  "!'Vr— 4«  *.ft  ■  ” 

7  ft  I  LOT  CCTTO'  O 

-  7  c  '  .  F.  *  9  <'■ 

- :  0  0  ?  c  0 . : 

-•  >  "ft  r;ft  c  t  i_  r 


TA8LE  9  -  SAMPLE  RUN  STREAM  FOR  CASE  B  (CONTINUED) 


DRAFT 


J.'  =  CA  ;D  »C  '*1 7  Z  COO  .  T20  0 

CM  »  s  r  .  •  •  4  mt  t*  *  *  *  *  »  ♦  *  « 

5'irjr?T*TAP?ii  . 

st,t  a?T!  '♦*"  r. 

l^T4C~,T£^ri.L->'l-,Y''tT'>  T”  = 

irTiCr»T3rf:?iLA^-yo;T»' 

i t r a c « r .*, r ’ •  t  ' ^  — v c  - r ~  .cv-'>.  t~  = 
*  IT  AC-  »  -  t  *0"  1  ~  ■  £T~°  ™’!  :»V'- 
A'TAC-,HYni'iT.Dr *  - r  u~ ,  :'•= 

Lj«D  »  -Y"I  '.T' 

:  n  as. 

si.i'DTi'ir ). 

C  '  -  Y  (  T  '  Pr  10,  ' 

P  r  «  T  ‘ 1  P  <  T  4  °  r  1 1  5  « 

CC  CY ( T  A P  £  1 1 t  CUT” U  T  )  * 


"  -  Wl’  i.  f  TAP£A  )  . 
r,PY!  T  APr  P  ,  CU  T="J  T  1  . 

-v.  r.  -,P  c  I  £  £ 

0 '••■  A  T  r  T"3T  c  r>r,r 

*SV;  ATHO  Tf’OF: 

x 0 '  i :  =  5.  *  •  . :  , 

—  1»7hF+3c»“1«0  3'.  *  l<  5  •  ■ 
1  .2  A">0C  , 

t>  ~  c  , -  T  •  - f)  £  *  n  ~  « 
&  •  A  f 

-1.  G  A  E”  *  c  ~  .-=>.7cr*''3  i 
1  •  0  AF+-0  A  , 

-1.  a:£+OA,-i  .lAr-  *cc  ♦ 

-1.77E  +  'jC*-l.<‘2r  ■*•?*!  » 
1.77=>0f  , 

-2 .OOT+O  f.  •  -!  .&lr  +  ^  ♦ 
C.O'TEl-fOf  ♦ 
ayfu:  =  a*'.'], 

-3.  = 

’.OIF*!!*  . 

—  1 . 0  f  *  0  f  ,  —  1  »  u  A 
l .  r  •:  =>  o  3 , 

—  p  •  n  'i  ~  *  a  '- , — 2  •  c .a 

.  o  A  r  >  r  c  , 

-3. 2  a  , *0&  . 

3.2  3- 

-2  a  « - : .  7°  r  ♦ :  a  « 

-2.3  ’p  +  o  a  . .si  7  ♦Of  » 


r  +  0  3  . 

“  ♦  0  A  , 

•■*2,1 


r  r 

3.  I  OPT  "(15=0.  TroT,f2):l*r,Tf^:T,n''Tr^rT. 

7r*r:4tl.  '  7  r  *  ^  r  • 

-.•3  nr,  -i  .  3--fc>  Cc.  1.  r*0r.2.  »Ar*n:,;.:or*^, 

if'  .3  ’  r*'e  »  -  •  A  •  ■+^  =  f'  '^  *  1  +  "  "  * 

1,1  -rr.rlt.  .nr*0r,.i.i  jr.QC,,  .l7r*Of  ,1  .c  ir  C  . 

7. «  »r* 2,_-  .i-  -.os . p.  i  ■  r+o  =  : ir*r  r*T.  r*rc . 

•  !. s.-  .073* "3,  --.o  ■»"♦?=.  .;*f'*rf  r 

.  t.  ,2rfH  ,_i  .7  ->r*r  6 . 1.3 -P  +  0  f-  ,l.?°F-0  A  •?.  ’A"*'  f  . 

■9.  T !*>!)*»  -1  f  ,  1.7'  r*0f  .3.51  ►  *0  c  .O.clr  +  r--» 


TABLE  10  -  SAMPLE  RUN  STREAM  FOR  CASE  C 


53 


'  K  P  t 1  2  —  1  *  ’  • 


i  7  ,■»  ' 


.  n  =>0*  •*.?  r>r 


-  ^  ^ 


-*  n  *7  ~  *  *  ' 
"«_•//  * 


_  t  i  i  r 


■  J  ■'  i 


-<r..7f.~+rJr-  ♦ 

I  .l*r«-0  7 ,1  7.1  .r :4'*n,-»  .?  A  7  -  0  s 

-  1.  1  J  7  , 

1 .5!?"*r*.l.ft'>r*«>'  ‘  f  "*'!  7«  1  .'j*  7 

-1. rcr+C 7, 

i. 7  nE  +  07. 1,7 7.  i  «-♦'>  ■» 
-1. 7  Cr*" 7  . 

a,\S'J'J  = 


*.l  ,7  I-*’  •". 

-5.  ur*c--, 

?.4C  r+07  1  4,>  ,-  7.^  2  .7  c""1-'1  / 

-2.  *  >>0  7  , 

&*ftftr*07*&«r;  I  -  +  ’,"  +  r  7#  ?  *  7!  *  j  7 

-4,  c;r+r.7, 

?.07E+O7,7.F7E*:7.7  .  '  .ft  ?"*r>7 

-7.Mr*D7. 

7.4cr+G’  .7. 4-7>C7 ,7.'-  =5=>0  7.5  .  c  ■  7 
-7.4-r*a 7. 

-:.f  "'■♦07,fe.f  r-r+C  7,G.f  w.n?,  s  .2  V  *0  7 


-  /t  .  ■  <  f  «.  •  <■ 


-  7 .  <> .-  f  +  ; 
-i.  '*e. ) 


1 


.  - 

-5  •  5  ;>t  ♦C  7 
-  7 . 7  7  :  7 


-  c  .  r  -. "  *  '.•  7 


7  L  +  ? 


'  T  . 


-■'  .Zt'+Z  ft. 

-1.04^*07, 
*0  7  . 

- :  *  t  *  r  *  o  7  • 


.  7  .  'i  rt  r  *.  - 


-'.?ir*07. 
--.7  7  '  ♦  C  7  . 
-'7.^kr  *^7. 
.ft~E *0  7i 

.  g^e*: 7. 


-ft.  ft"r«-C 


T  *A  y  rl2f?  .  »=,TCH?C=1  .  7  7  .  'TCHcr-1,'7, 

p?pMco  =  tf  C..  -  “PT -1  .  nT:'.^  *■’  ”  T  *•’ T  ~i  ♦ 

~  t>TL7)7  A CT  * r  '"r 


1ft  1  7  .  ~  '- 
:  ft  1 7 .  ^  ■’ 

i  nc  o .r>o r 

1302CC.C 

V  *. T H  TEST  E "D 

T 5V  A  TM?  *:c>j-  r-a.  n 

r  =  i  let  ACTirsfis 


-  7  *  .  n 


ion :  .r  +  "  r 
•  l  o  o  o p  t> .  n 

■.•  ■  t n  nc  r  t  i  — 

„  *■  *—  — 


TABLE  10  -  SAMPLE  RUN  STREAM  FOR  CASE  C  (CONTINUED) 


J'fJCArD«Cv  !  "T  ~C  :■  *  T2''  j  »r4 . 

Ch*A'4';T*Ta  ‘•p.  ****-..***■■.. 

’•TJflST.T^r,  l  . 

'-::,jroT,raf 

’-:faUC?T*Ttrr,:  .*P  r. 

\  T  T  A  C  “  .  T  4  c  c.  1 .  L  *  K  -.  y  n  A  T  -•  .  c  v  Z  -  ,  I r  = 

aTTAC--'»TlPFp,!.;u-~iY-.-rr.“v  =  -j,  i^  = 

A  AT  AC  .,TJD  F  *  ,!_*.  ✓  i  v~i  -rr  .'‘Yr’',  T  "  = 

aTTAC'"»=r  A-* c  •AT-o  !l'"",T  ''r 
att  Acu*HYrif-T,  •->:*••  i^Yl &=»!“=  ' 

L  iTf-'iDr:n 

: :  \  *  r  . 

5-:i,'irjT53r.! -• ) . 

C  ’PY  (TAT"]  .  :‘i'T~UT)  . 

^vir^fTspcin. 

C  '°Y(  Tier  1  1  ,-'uT  =  ',Tl  » 

C'PY  <  TA^o  .CUT3 ‘J T  )  . 
r»>r  o=  F  I|_" 

~  A  TH  TEST  ~  ~-~,r  C*"r  ' 

'SWATH*'  Tn'ic  =  r.  (  1'  =T  «  (  2  '  =t  •  T  <'•'  >0  ' 

AXPUT=  P*:.T  , 

-i.  - p i-i  .:Pr*r,p .  - 1. ’  r*-  ' -  r r  t  — +  t  ■-  .i. :  ■"p*r:  .  r  - r*rc 

1.5^P+0~  , 


-p.4P.T*P5,-4.GOr^O=.-3.5«-E>T5 

-1 .  3  4T+:  c  ,  r*2~  , .  2  ?<->  ~  E 

l.QA;>0f  • 

-1.5iE+Pb  +bh*-7.r. •■•  -  +  CC 

1  . 

-!.77r+,’n,-l  .  .i  1  £  *  0  .-0.4  :t>r.* 

1  ,7Tr*rc- , 

-i.  p  f .ir,r‘np  - 

A  YRIJJr  4*C.'« 

-4.  ?  =  ."♦;■  T  ;  r 

4.23E+CE, 

-1.3?E>P6«-!  •  ?  I  E  *T  e  « -  i  .  .'E+  j  c 
!  .7  7  ->3  a  , 

-Z.  A  7E>r  P  .  nl  -"*CP.  •  —  I  •  7  F*  : 
..0  7  3  ♦  C  , 

-7,c-'G  +  n;  ,-~.2  T*  Tf. 

z .a o r*-f. . 

.  1  71>G  t  .-1  .*K  r*cr  ■-»:>  "  ^ 

”'.17r-*-jp« 

-1.3~'**oc  .-]  .Zl” +2  c-  i. ^  ' 

U?---*C0  * 


=  .  3 . 1  '  r  ♦  C  ; 


7  .  1  -  ‘  ♦  P  r 


?.-r  *':t.  •  4.2cr  +  'l  = 


♦  G 


pr-  +  n~,  T.  4  •-  r  *  !7 r: 


•r*Pc.c.?''r'  +  G: 


7  r  C  .  ’  .  K  '  ♦ : 


7 « "♦'? c , 7t  p 


1  '  C  ♦  C  A  .  ;  .  i  --  S  ♦  :  - 


?P  -  V  ,<=.  ,o.3fi  -4-0  = 


f  .  ?  r  L  * r  C  »  r 


.  r  o  r  4.P  =  .  i 


.  4 1  e  ■*.  - c , : 


7.2Pr  +  r'~.7 


.  «■  ,'J  F  *  p.  c  7 


} .  i  •.'*7  o 


1.7 


2  .  ?  7  r  ♦  0  c  .? 


!  .  P  4  E  *  •'  f  .  1 


i  .t:r+o a 


7'f  ♦’  - 


•]>>'  A 


1  "  *  r  P 


u  o  f  *  r  ^ 


TABLE  11  -  SAMPLE  RUN  STREAM  FOR  CASE  D 
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O  lT* 


.  c  c  -  *  ■ 


'j:.=  2  » 

r  .  ?r  «  c.  1  ??+  •, -«4.=  . 

GT*0  E  « 

r+ c  c. .  r .  a  4  r  ♦  r  c  .  ?  .  r-  r  f 2  f ,  l .  7  i  '  «• z  ~ 
-2.  ?:F*0‘.t 

■s  .7."  *r  +  2^ .0 . *  ir  *c  <i  .g <'r  +  : ' 

. 3CF*0  £  . 

"TOC  . 

2 .7?f+:g  *7.’3ftF* :  c ,? .  i  7r*  •  -  .c :  i r  *r  - 

-a.7-r+3c  , 

-G. 1 3E-1 c. 
ir.?!l~r  7  *  C  .  C  , 

s.9l>rfi  7E*>'-.5  .2~F-CGt  7  .??"♦!!£ 

-  -  •  °  1  r  ♦ c » 

7 


1  7  ^Z'7 

c  lz* 7  7  *  3  .  ’)  3  •”  *  ft  7 


-4.  "c>: 7 . 

-4. 77E+: 7  . 

,"  =  -*■  Z  7  •?  ,Q  ^■T‘-r' 7  •  2  .c?r  '■Z'7 

-*.b*l*C  7. 

4  .?7-  +  0  7»’,.F*;r*"7»7.,rr*.17.'>.-"ir+--> 

-4.  2  7r*0  7  . 

t'*4v=i21  .,-TC-c  -=1  »  7  7  t  ■'  T  Cr!  "  !'c  1  ,  "  t  , 

3^?'-<EO  =  lc'>..  oq=-*PO  =  lB0.  •  DT=7.r*  ■’2'^T»'T  = 
-  ot|_!}t  :CT  r"’  ? 

P  1 7 . c  q  7  -7-.0 

ic  17. 71  r  “  -  7  r  .  0 

_  ■»  c*  r*  Ta.*'* 

•  •  o  •  L  ♦  -  ^ 


-i .  7ir+ 1 £ « -  ~ .  i  - r * o c 

-c.?rroc  «-'.'i~r+0c.-Q.0''r*rr 

-c.c  77 +  Z.'  .-7.!  7r  +  0  «■ ,  -7.=  cr+r  ‘ 
-7.-cE  +  0-  .-4."t2r*0G.-5.crr*Pc. 

-7.'7:r  +  r;F,--.22!r*Cc.-5.7',r+Cc 

-1.P7C  ,«tc.P,.?,7'r^  7 

- 1 .  ?  rr  «•  ■■ ■  ’ .  -  & .  i  «=•  f  *o  7 .  -4 .  <+ r- 1  ♦ r- 7 
-?.°?E>n7,-7.  =  6r.kQ7,_4.o  =  F>r7 
-2.74E*::*-T.7cr+C;7.-T.I',r*"' 


$* c t h  test  r*r 

•»>y  4TM[)  '  5|iMr-*  •  IT’  S 

”  °TL7T  •'•  CT  T  " ’  S 
-1  SC  0  .r+'7  ~ 

-icc  :r-s • 

re.'  ** e  c  t  i  r 


TABLE  12-DEFINITION  OF  FORTRAN  SYMBOLS  APPEARING  IN  NAMELIST  /SWATHD/ 


r 


DS 


FORTRAN 

Symbol 

Description 

Type 

Units 

AE 

Absolute  error  criterion  for  range-KUTMER 

Real 

- 

AKRUD 

integration 

K  rudder  moment  array 

Real 

nt-m 

ANRUD 

N  rudder  moment  array 

Real 

nt-m 

ARAFIN 

Areas  of  fins 

Real 

m2 

ARATIO 

Rudder  sizing  ratio,  to  standard  rudder 

Real 

- 

AXRUO 

X  rudder  force  array 

Real 

nts 

AYPUD 

Y  rudder  force  array 

Real 

nts 

BRDARM 

Aft  rudder  pitch  moment  arm 

Real 

m 

BRSPAN 

Maximum  bow  rudder  extended  length 

Real 

m 

BRUDPO 

BRUDSO 

CLFIN 

Port  bow  rudder  position,  initial  time 
Starboard  bow  rudder  position,  initial  time 
Lift  coefficient  of  fins 

Real 

/rad 

DT 

Step  size  for  integration  procedure 

Real 

s 

DTMIN 

Minimum  step  size  for  KUTMER 

Real 

s 

EPS 

Relative  error  criterion  for  KUTMER 

Real 

- 

FINANO 

integration 

Fin  angles  at  initial  time 

Real 

rad 

HULPIK 

Difference  in  draft  from  design  condition. 

Real 

m 

I  DOF 

deeper  draft=positive 

Degrees  of  freedom 

Integer 

IOPTN 

Integration  technique  and  optional  pages  3-5 

Integer 

? 

IXX 

Ship  moments  of  inertia 

Real 

kg-nu 

IXY 

Ship  products  of  inertia 

Real 

kg-nu 

IXZ 

Ship  products  of  inertia 

Real 

kg-nu 

IYY 

Ship  moments  of  inertia 

Real 

kg -itu 

IYZ 

Ship  products  of  inertia 

Real 

kg-nu 

III 

Ship  moments  of  inertia 

Real 

kg-nu 

KPDOT 

Dimensional  coefficient  K* 

D 

Real 

kg -nr 

KRDOT 

r 

Dimensional  coefficient  K* 

r 

Real 

2 

kg-m 

KVDOT 

Dimensional  coefficient 

Real 

kg 

LENGTH 

Ship  length 

Real 

m 

MASS 

Ship  mass 

Real 

kg  2 

MQDOT 

Dimensional  coefficient  M* 

Real 

kg-m 

MSUBTH 

Restoring  derivative  MQ  due  to  buoyancy 

Real 

nt-m/rad 

MSUBZ 

Restoring  derivative  M2  due  to  buoyancy 

Real 

nt-m/m 

MWDOT 

Dimensional  coefficient  M- 

w 

Real 

kg 

NPDOT 

Dimensional  coefficient  N* 

n 

Real 

kg-m2 

NPRINT 

r 

Number  of  DT  increments  in  print  loop 

Integer 

- 

NRDOT 

Dimensional  coefficient  N- 

r 

Real 

2 

kg-m 
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TABLE  12  (continued) 
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FORTRAN 

Symbol 

Description 

Type 

Units 

NRUN 

Simulation  run  number 

Integer 

- 

NUMRAL 

Number  of  data  grid  rudder  angle  points 

Integer 

NUMUUV 

Number  of  data  grid  velocity  points 

Integer 

- 

NVDOT 

Dimensional  coefficient  N^ 

Real 

kg 

PDIA 

Prop  diameter 

Real 

m 

PHIO 

Initial  0  value 

Real 

- 

PNT 

Optional -debugging  output,  9  =  hydro  array 

10  =  intermediate  force-moment  time  histories 

Logical 

- 

PO 

Initial  P  value 

Real 

- 

PROARM 

Prop  X,  Y,  Z  coordinates 

Real 

m 

PRPMPO 

Initial  port  propeller  rpm's 

Real 

rpm 

PRPMSO 

Initial  starboard  propeller  rpm's 

Real 

rpm 

PSIO 

Initial  y  value 

Real 

- 

PTCHPO 

P/D  ratio  -  Dort  at  T=0 

Real 

- 

PTCHSO 

P/D  ratio  -  starboard  at  T=0 

Real 

- 

qo 

Initial  Q  value 

Real 

- 

RALFA 

Grid  rudder  angle  array  for  rudder  forces 

Peal 

deg 

RO 

Initial  R  value 

Real 

- 

RUDANO 

Rudder  angles  at  initial  time 

Real 

rad 

RUDARM 

Aft  rudder  pitch  moment  arm 

Real 

m 

SCFACT 

Ship-Model  scale  factor 

Real 

- 

THETAO 

TITLE 

Initial  0  value 

Title  of  run 

Real 

TMAX 

Maximum  time  for  run 

Real 

s 

TO 

Initial  time 

Real 

s 

UO 

Initial  U  value 

Real 

- 

UUVEL 

Grid  velocity  array  for  rudder  forces 

Real 

m/s 

VO 

Initial  V  value 

Real 

- 

XC 

X  coordinate  of  center  of  gravity 

Real 

m 

XFIN 

X  coordinate  of  fin  center  of  pressure 

1  =  bow  starboard,  2  =  bow  port,  3  =  aft 
starboard,  4  =  aft  port 

Real 

m 

XRUD 

X  coordinate  of  aft  rudder  center  of  pressure 

Real 

m 

XUDOT 

Dimensional  coefficient  X^ 

Real 

kg 

YC 

Y  coordinates  of  center  of  qravity 

Real 

m 

YFIN 

Y  coordinates  of  fin  center  of 
pressure 

Real 

m 

.  2 
kg-m 

YPDOT 

Dimensional  coefficient  Y- 

D 

Real 

YRDOT 

r 

Dimensional  coefficient  Y< 

r 

Real 

kg-m^ 

YRUO 

Y  coordinates  of  aft  rudder  center  of 

Dressure 

Real 

m 
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TABLE  12  (continued) 


FORTRAN 

Type 

Symbol 

Description 

Units  j 

Dimensional  coefficient  Y* 

Real 

kg 

p  ''  I 

Z  coordinates  of  center  of  gravity 

Real 

m 

Z  coordinates  of  fin  center 

of  pressure 

Real 

m 

i  1  =  bow  starboard,  2  =  bow  port,  3  =  aft 

starboard,  4  =  aft  port 

Real 

.  2 

ZQDOT 

Dimensional  coefficient  Z^ 

kg-m 

ZRUO 

Z  coordinates  of  aft  rudder 

center  of 

Real 

m 

ZSUBTH 

pressure 

Restoring  derivative  ZQ  due 

to  buoyancy 

Real 

nts/rad 

ZSUBZ 

Restoring  derivative  lz  due 

to  buoyancy 

Real 

nts/m 

ZWOOT 

Dimensional  coefficient  Z^ 

Real 

kg 
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the  Namelist  input  data.  The  FORTRAN  variable  NRUN  should  be  assigned  a 
unique  value  for  all  runs  to  be  submitted  at  one  time.  Only  ID0F=5,  and 
ID0F=6  modes  of  operation  have  been  tested.  During  the  course  of  this  study, 
several  values  of  DT  were  used.  A  value  of  3.0  results  in  a  very  few  inter¬ 
val  halving  operations  in  KUTMER,  and,  hence,  can  be  considered  optimum.  A 
smaller  value  of  DT  may  be  desirable  due  to  print  out  considerations.  Current 
ly,  state  variable  time  histories  are  printed  at  integer  multiples  of  DT, 
every  NPRINT*DT  seconds.  If  a  sampling  interval  of  3.0  seconds  is  undesirable 
DT  should  be  changed. 

Pilot  Action  Data 

The  control  system  allows  the  programmer  to  read  in  up  to  ten  control 
actions  via  pilot  action  data  cards.  The  system  is  designed  so  that  the  pro¬ 
grammer  can  choose  from  twelve  controlled  items  and  sixteen  state  controlling 
variables.  The. system  allows  for  both  single  instance  or  constant  monitoring 
of  the  controlling  variables.  For  example,  assume  the  aft  rudder  is  tc  te 
turned  to  -20°  if  <  10°.  If  single  instance  is  chosen,  this  action  will  te 
performed  the  first  time  ip  <  10° ,  Linder  constant  monitoring,  however,  ? 
action  will  be  performed  every  time  <  10°. 

In  the  present  system,  the  controlled  item  must  have  a  numerical  end 
value  read  in.  There  is  no  provision  for  reducing  a  variable  by  ten  percent, 
for  example.  All  pilot  action  data  items  are  in  10  column  fields  right 
justified*.  The  first  card  of  data  indicates  how  many  pilot  action  cards 
are  to  be  read.  The  pilot  action  cards  have  the  following  format: 


♦Floating  point  numbers  need  not  be  right  justified  provided  the  decimal 


point  appears  explicitly. 
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J1 


CV1(I)  J2 


CV2  ( I ) 
31-40 


i'iw 


3 


Columns  1-10  11-20  21-30 


CV  3 ( I )  JACT(I) 

41-50  51-60 


where 

J1  =  index  of  state  vector  controlling  variable 

CV1  =  numerical  value  of  J1  when  control  is  to  be  activated 

J2  =  controlled  item 

CV2  =  final  value  of  controlled  item 

CV3  =  rate  of  change  of  controlled  item 

JACT  =  1  =  single  instance,  0  =  constant  monitoring. 

The  definitions  of  possible  values  for  J1  and  J2  are  given  in  Table  13. 

The  units  for  CV1,  CV2,  and  CV3  are  MKS  units.  Angles  are  in 
degrees.  If  a  control l-'ng  variable  condition  is  to  be  less  than  some 
CV1  value,  e.g.  If  YE  <-  100  m,  the  index  of  the  controlling  variable 
should  be  assigned  a  negative  value.  C0NTRL  is  called  during  each  call  to 
DAUX.  00NTRL  checks  to  see  whether  the  controlled  item  needs  to  be  changed, 
and  if  so  performs  the  changes  at  the  desired  rate.  When  the  control  item 
reaches  its  final  value,  C0NTRL  checks  JACT  and  if  JACT=1,  the  control  item 
is  turned  off.  Otherwise,  that  controlled  item  is  continuously  monitored 
for  possible  further  actions.  If  multiple  actions  have  been  programmed  for 
a  controlled  item,  an  error  message  will  be  printed  out  and  the  last  action 
performed. 

Rate  Dependent  Coefficients 

The  pilot  action  data  is  followed  by  rate  coefficients.  Currently, 
values  are  not  available.  The  values  -1000. E  +  00  and  -  100000.  indicate 
to  the  subroutine  INPUT  that  rate  coefficients  are  not  to  be  read. 
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TABLE  13  PILOT  ACTION  DATA  DESCRIPTIONS 


.T^T) 

m 


j-aj  a 


INDEX 

STATE  VECTOR 
CONTROLLING 
VARIABLE 

CONTROLLED  ITEM 

1 

u 

Starboard  propeller  RPM 

2 

V 

Port  propeller  RPM 

3 

w 

Starboard  propeller  P/D  ratio 

4 

p 

Port  propeller  P/D  ratio 

5 

q 

Forward  starboard  rudder  position 

6 

r 

Forward  port  rudder  position 

7 

Roll  angle 

Aft  starboard  rudder  angle 

8 

Pitch  angle 

Aft  port  rudder  angle 

9 

Yaw  angle 

Bow  starboard  fin  angle 

10 

XE 

Bow  port  fin  angle 

11 

YE 

Aft  starboard  fin  angle 

12 

ZE 

Aft  port  fin  angle 

13 

Distance 

travelled 

14 

Speed 

15 

Drift  angle 

16 

Time 
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DESCRIPTION  OF  OUTPUT 

An  example  of  printout  is  presented  in  Appendix  C.  A  brief  descrip¬ 
tion  of  the  output  follows.  The  title  card  and  all  data  read  via  Namelist  is 
printed  at  the  beginning  of  the  output.  This  output  appears  as  the  FORTRAN 
symbol  followed  by  the  assigned  value.  After  the  Namelist  is  printed  the 
requested  pilot  actions  are  listed.  The  hull  hydrodynamic  data  are  next 
printed  if  PNT(9)=.TRUE. . 

Following  the  above  output,  up  to  5  pages  of  time  histories  are 
printed.  These  contain  time  histories  of  the  state  variables  and  time  histo¬ 
ries  of  the  control  surface  deflections.  Pages  2-4  are  printed  if  I0PTN(2) 
is  set  to  1.  These  pages  consist  of  time  histories  associated  with  the  rud¬ 
der,  fin,  and  propulsion  subsystems.  The  printout  starting  with  the  Namelist 
and  ending  with  page  five  is  repeated  for  each  run. 

Intermediate  forces  and  moments  are  printed  following  the  paged  output 
for  all  runs  (if  PNT(10)=.TRUE. ) .  These  data  contain  a  time  history  of  each 
force  component  computed  in  DAUX.  A  complete  table  of  X,  Y,  Z,  K,  M,  N  forces 
and  moments  is  given  at  each  time  step.  This  option  is  useful  for  debugging 
purposes,  as  unrealistic  ship  motions  are  often  caused  by  improper  program¬ 
ming  of  a  specific  force  subroutine. 

A  time  history  of  time  steps  used  by  KUTMER  is  printed  at  the  end 
of  each  computer  job  submitted  if  PNT(10)  =  .TRUE. .  Optimum  integration  time 
steps  can  be  selected  based  on  an  examination  of  this  output. 

The  next  section  contains  descriptions  of  various  subroutines. 
Sufficient  detail  is  presented  so  that  new  force  and  moment  models  can  be 
easily  developed. 
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SUBROUTINE  DESCRIPTIONS 

This  section  contains  a  short  description  of  the  main  program  and 
the  associated  subroutines.  Each  description  contains  the  purpose  of  the 
subroutine,  how  to  call  it,  parameter  descriptions  and  clarifying  remarks.  Ex¬ 
cept  for  the  main  program,  the  subroutines  are  listed  in  alphabetical  order. 
Program  SWATH 

Purpose  -  Main  driver  for  simulation.  Calls  INPUT  to  initialize 

variables,  calls  KUTMER  in  integration  loop,  calls  OUTPUT, 

and  checks  to  see  if  all  runs  have  been  completed. 

Subroutines  called  -  INPUT,  KUTMER,  OUTPUT 

Remarks  -  1)  SWATH  calls  INPUT  to  initialize  variables 

2)  SWATH  calls  KUTMER  in  a  do-loop  of  order  NPRINT,  thus 
obtaining  output  values  in  increments  of  DT  x  NPRINT, 
e.g.,  DT  =  1,  NPRINT  =  3,  output  is  in  3  sec  intervals 

3)  SWATH  automatically  terminates  the  run  when  time  >  TMAX 

4)  Each  run  must  have  a  different  NRUN  value,  a  repeated 
value  or  a  negative  value  stops  the  program  and  starts 
intermediate  printout. 

5)  SWATH  calls  OUTPUT  in  the  integration  loop  to  load 
values  into  output  arrays.  OUTPUT  is  also  called  outside 
the  loop  at  the  end  of  a  run  to  print  out  the  last  page 
of  output  that  will  not  necessarily  be  a  full  forty 
lines. 

6)  Output  is  printed  one  page  (40  line)  at  a  time.  The 
data  is  loaded  into  page  arrays 
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Subroutine  ADDFOR 

Purpose.  To  calculate  the  hull  forces  due  to  first  and  second  order 
rate  dependent,  terms 
Usage.  CALL  ADDFOR 

Remarks.  1)  The  first  order  equation  used  is 

ADDFM(J)  =  RTMRX(J,K)  *  Variable  (K)  •  dimensional izing 
factor 

2)  The  second  order  equation  used  is 

ADDFM(J)  =  ADMTRX(J,K,L)  *  Variable  (K)  •  Variable  (L) 

■  dimensional izing  factor 

3)  RTMTDX  and  ADMTRX  are  the  derivatives  of  the  forces 
and  moments  with  respect  to  the  volocities,  e.g.  Ny,  Xu 

V‘ 

4)  The  dimensional izing  factors  are: 

*sp(ship  length)2  for  J=l,2,3  (X,Y,Z) 

%p( length) 3  for  J=4,5,6  (K,M,N) 

Subroutine  AIRFOR 

Purpose.  To  calculate  the  aerodynamic  forces  and  moments 

Usage.  Call  AIRFOR 

Remarks.  l)  Estimates  show  that  the  aerodynamic  forces  can  be 
neglected  and  all  forces  and  moments  are  set  equal 


zero. 


Subroutine  BOYFOR 

Purpose,  to  calculate  forces  and  moments  due  to  bouyancy  and  gravitation 
Usage.  CALL  BOYFOR 

Remarks.  1)  Equations  used  for  buoyancy  are 

Z  -  zhh  +  v 

K  ■  ■  tan  0  *  Z2 

M  =  Mhh  +  M06 

2)  Zh,Z0,Mh,Me  are  assumed  constant  for  small  motions  and  given 
in  the  body  coordinate  system 

3)  The  gravitational  moments  are  calculated  as  weight  *  moment 
arm,  where  the  moment  arm  is  dependent  on  the  distance  Xg.Yg, 
or  and  the  ship  orientation. 


Program  BLOCKD 

Purpose.  Initializes  all  constants  and  ship  state  variables  appearing  in 
named  common  blocks 

Remarks .  1)  All  variables  have  been  initialized  to  perform  straight 

ahead  flight  at  20  knots. 
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Subroutine  CONTRL 


Purpose,  to  determine  whether  a  pilot  action  is  necessary  and  if  so,  to 
perform  the  operation  at  a  specified  rate. 

Usage.  CALL  CONTRL  (T.Y) 

Description  of  Parameters. 

T  =  time 

Y  =  state  vector,  see  INPUT  for  description 
Remarks .  1)  'control  actions  are  read  in  on  cards.  See  users  guide  section 

description  of  parameters  and  Table  13. 

2)  Controlled  items  are  fin  and  rudder  angles,  bow  rudder  position 
and  propeller  pitch  and  rpm. 
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Subroutine  DAUX 

Purpose,  to  compute  the  right  hand  side  of  the  differential  equations 
of  motion.  (See  Figure  11  for  Flow  Chart). 

Usage.  CALL  0AUX(T,Y,F) 

Description  of  Parameters. 

T  -  time 

Y(l)  -  U,  velocity  along  body  x  axis  (m/sec) 

Y(2)  -  V,  velocity  along  body  y  axis  (m/sec) 

Y ( 3 )  -  W,  velocity  along  body  z  axis  (m/sec) 

Y(4)  -  P,  angular  velocity  about  body  x  axis  (rad/sec) 

Y(5)  -  Q,  angular  velocity  about  body  y  axis  (rad/sec) 

Y(6)  -  R,  angular  velocity  about  body  z  axis  (rad/sec) 

Y(7)  -  <j>,  roll  angle  (rad) 

Y(8)  -  9,  pitch  angle  ( rad) 

Y(9)  -  ip,  yaw  angle  (rad) 

Y(10)  -  X^,  x  position  in  fixed  coordinate  systems  (m) 

Y ( 1 1 )  -  Yg,  y  position  in  fixed  coordinate  systems  (m) 

Y ( 12 )  -  z  position  in  fixed  coordinate  systems  (m) 

Y ( 13 )  -  S,  distance  traveled  (m) 

F( 1- 13)  -  derivatives  corresponding  to  the  above  variables 
Subroutines  Called.  CONTRL,  FINFOR,  RUDFOR,  HULFOR,  AIRFOR,  BOYFOR, 
PROFOR,  ADDFOR,  LINQ 

Remarks.  1)  The  6x6  mass/inertia  matrix  is  computed  from  the  ship 
mass,  moments  and  products  of  inertia,  and  added  mass- 
inertia  terms. 
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FIGURE  11  FLOWCHART  FOR  SUBROUTINE  DAUX 
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2)  The  force  vector  is  determined  by 

A  -  calling  the  separate  force  subroutines 
B  -  calculating  the  terms  due  to  space-body  transformation 
C  -  calculating  terms  due  to  center  of  gravity  not  origin 
D  -  summing  all  forces 

3)  DAUX  computes  the  space-body  transformation  matrix  based  on 
the  ship  orientation.  This  matrix  is  used  to  calculate  the 
fixed  space  acceleration  terms  from  the  body  axis  acceleration 
terms  returned  from  LINQ. 

4)  If  degrees  of  freedom  are  less  than  6,  DAUX  uses  a  constraint 
equation  to  evaluate  the  appropriate  terms. 

Function  DRAAGG 

Purpose,  calculates  hull  drag  at  a  given  speed 

Usage.  DRAAGG(V) 

Description  of  Parameters. 

V  -  Speed  (m/s) 

Remarks.  1)  Drag  is  calculated  as  EHP  x  1450.1 

V(m/s) 

2)  EHP  is  computed  by  interpolating  a  data  base  of  EHP  values 
for  speeds  from  0  to  28  knots  with  grid  points  at  integer 
values  of  speed. 

Subroutine  EULER 

Purpose,  performs  EULER  integration  on  state  vector  Y 

Usage.  CALL  EULER  (T,DT,Y,NSTATE) 

Subroutines  Called.  DAUX 
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Description  of  Parameters. 

T  -  time 

DT  -  step  size 

Y  -  state  vector  (See  INPUT) 

NSTATE  -  number  of  variables  in  Y 

Remarks.  1)  DAUX  is  called  to  obtain  the  time  derivatives  of  state  vector 
Y 

2)  Y  =  Y  ,  .  +  DT  *  F  where  F  is  the  time  derivative  vector 
'  new  old 

3)  EULER  has  no  error  check  and  is  a  purely  mechanical  subroutine. 
If  too  large  a  DT  is  used,  there  will  be  major  errors  in  the 
integration  after  many  time  steps. 

Subroutine  FINFOR 

Purpose,  to  calculate  the  total  fin  forces  and  moments 
Usage.  CALL  FINFOR 
Subroutines  Called.  TOTVEL 

Remarks.  1)  TOTVEL  is  used  to  calculate  the  total  inflow  velocity  at 
each  fin 

1  2 

2)  Fin  left  is  calculated  theoretically  as  Lift  =  ^pAU  C[_ 


where  p  =  density  of  water 
A  =  fin  chord  •  span 
U  =  total  inflow  velocity  Vu^+ 

=  lift  coefficient 

3)  Fin  drag  approximated  as  Lift/7 

4)  Center  of  pressure  assumed  at  quarter  chord  and  half  span. 

5)  Effects  of  side  flow  (v)  have  been  ignored 


Subroutine  HULFOR 


Purpose,  to  calculate  the  hydrodynamic  Y,K,N  hull  forces  as  a  function 
of  speed,  roll  angle,  yaw  rate,  and  drift  angle. 

Usage.  CALL  HULFOR 

Subroutines  called.  DRAAGG,  HYDINT 


Remarks.  1) 

Equations  used. 

at  a  given  speed 

Y  -  V  +  Yr,r' 

+  YJ 

K  -  KbS  +  Kp , r 1 

+  KJ 

<p 

N  =  Ng6  +  Nr,r' 

+  N  <j> 

2) 

HYDINT  and  HYDPAR  perform  a  four-dimensional  interpolation 

to  obtain  the  Y 

,K,N  forces  as  a  function  of  speed,  drift  angle, 

yaw  rate,  and  roll  angle. 

3) 

DRAAGG  calculates  the  X  force 

Subroutine  HYDINT 

Purpose,  performs  a  4  dimensional  interpolation  for  hull  hydrodynamic 
forces  based  on  data  curves.  Independent  variables  are  speed,  drift  angle, 
yaw  rate,  and  roll  angle. 

Usage.  CALL  HYDINT 
Subroutine  Called.  HYDPAR 

Remarks.  1)  Hydrodynamic  values  must  be  stored  in  array  HY(4,3213)  in 
a  linear  fashion.  The  first  index  represents  the  X,Y,K,N 
data.  Array  XG(4,20)  contains  the  data  base  grid  values. 

The  first  index  represents  the  four  independent  variables: 

1  *  roll  angle,  2  =  drift  angle,  3  =  yaw  rate,  4  =  speed. 
Array  LG(4)  contains  the  number  of  base  grid  values  for 
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each  parameter.  The  indexing  of  HY ( Y , 32 13 )  is  as  follows: 
for  Array  HY(I,J),  the  Jth  value  corresponds  to  the  I  force 
at  the  location 

J  =  IR  +  LG (1 ) F  ID-1  +  LG(2)  {IY-1  +  LG(3) ( IS-1) }] 
where  IR,  ID,  IY,  and  IS  represent  the  grid  point  index  in  the 
XG  array  for  roll,  drift,  yaw  and  speed,  respectively  and  LG(k) 
equals  the  number  of  grid  points  for  that  parameter. 

Subroutine  HYDPAR 

Purpose,  to  recompute  the  interpolation  coefficients  for  HYDINT  if  the 
data  cell  is  not  entered  when  HYDINT  is  called 

Usage.  CALL  HYDPAR 

Remarks .  1)  HYDINT  will  call  HYDPAR  automatically  if  necessary. 

2)  The  coefficients  calculated  by  HYDPAR  are  based  on  the  cell  grid 
points  and  the  data  base  grid  points  need  not  be  uniformly  spaced. 

Subroutine  HYRITE 

Purpose,  optional  listing  of  hydrudynamic  hull  coefficients  in  tabular 


Usage.  CALL  HYRITE 

Remarks.  1)  Called  from  INPUT  if  PNT(9)  =  .TRUE. 

Subroutine  INPUT 

Purpose.  Reads  in  NAMELIST  data,  hydrodynamic  hull  data,  pilot  control 
.1.4  and  rate  dependent  terms  and  converts  initial  conditions  to  MKS  units. 
^  '.all  INPUT(Y ,  TIMEE) 

■  r  Parameters. 

■  -  Snip  state  vector 

.  ’  velocity  along  body  x  axis  (m/sec) 


73 


2  =  V,  velocity  along  body  y  axis  (m/sec) 

3  =  W,  velocity  along  body  z  axis  (m/sec) 

4  5  P,  angular  velocity  about  body  x  axis  (rad/sec) 

5  =  Q,  angular  velcoity  about  body  y  axis  (rad/sec) 

6  =  R,  angular  velocity  about  body  z  axis  (rad/sec) 

7  =  4),  roll  angle  (rad) 

8=0,  pitch  angle  (rad) 

9  =  4),  yaw  angle  (rad) 

10  =  XE,  X  coordinate  in  fixed  space 

11  =  YE,  Y  coordinate  in  fixed  space 

12  =  2E,  Z  coordinate  in  fixed  space 

13  =  S  -  distance  travelled 

TIMEE  =  time 

Remarks .  1)  All  initial  conditions  and  values  that  may  change  due  to 

configuration  changes  can  be  read  in  via  NAMELIST. 

2)  Hull  data  read  in  from  tapes.  These  must  be  attached  in 
the  following  manner.  Tape  1  =  Y  data,  Tape  2  -  K  data, 
Tape  3  =  N  data. 

3)  Pilot  action  data  is  described  in  the  User's  Guide  section 
and  Table  13. 

4)  Rate  dependent  terms  are  read  in  as  follows 
Example  -  for  first  order  terms 

Ny  =  RTMTRX(6,2) ,  Kp  =  RTMTRX(4,4) • ■ • 1  =  X,U  2  =  Y,V 
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.  / '  < 

SJ 

3  =  Z,W  4  =  K,p  5  =  M,q  6  =  N,r 
for  second  order  terms 
Nvr  =  ADMTRX  (6,2,6)  and  (6,6,2)  as  in  the  above 
correspondence 

Subroutine  KUTMER 

Purpose.  Performs  Runge-Kutta-Merson  integration. 

Usage.  CALL  KUTMER  (ND,  T,  H,  VO,  EPSE,  A,  HCX,  FIRST) 

Subroutines  Called.  DAUX 
Description  of  Parameters. 

ND  -  number  of  variables  in  YO 

T  -  time 

H  -  step  size 

YO  -  state  vector,  enter  with  initial  values  at  T,  returns 

values  at  T  +  H 

EPSE  -  relative  error  criterion  for  YO  components 

A  -  absolute  error  criterion  for  YO  components 

HCX  -  smallest  step  size  to  be  used  by  KUTMER 
FIRST  =  0.  when  KUTMER  is  entered  for  first  time  or  with  changed 
step  size 

=  1.  when  entered  with  same  H  to  continue  integration  in  a 
sequence 

=  2.  When  error  criterion  cannot  be  met  and  H  <  HCX 
Remarks.  1)  KUTMER  calls  DAUX  five  times  to  predict  endpoint  value  of  Y 
2)  If  error  criteria  are  not  met  KUTMER  automatically  halves  the 
interval  and  repeats.  If  halved  interval  becomes  less  than 
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H/32,  KUTMER  stops  and  prints  out  error  message 


Subroutine  LINQ 

Purpose.  Solves  for  X  where  AX  =  B 
Usage,  CALL  LINQ  (A,B,N,KS) 

Description  of  Parameters. 

A,B  =  input  matrix,  LINQ  solves  for  X  and  replaces  B  by  X 
N  -  number  of  equations,  dimension  of  B 
KS  -  error  code,  1  =  singular  matrix,  0  =  normal 
Remarks .  1)  X  represents  the  acceleration  matrix  used  by  KUTMER 

2)  A  =  mass-inertia  matrix  obtained  from  DAUX 

3)  B  *  force  matrix  obtained  from  DAUX 
Subroutine  NGHBR1 

Purpose.  To  find  the  nearer  neighbor  less  than  or  equal  to  the  desired 
value  in  the  Array  V 

Usage.  CALL  NGHBR1  (V,  ARG,  II,  N,  IER,  RATIO) 

Description  of  parameters. 

V  -  variable  array 

ARG  -  search  argument 

II  -  Location  such  that  V(II).  LE.  ARG.  LE.V(II+1) 

N  -  number  of  points  in  V 

IER  -  0  in  range,  otherwise  out  of  range 

RATIO  =  (ARG  -  V(II ) )/  (V(II+1)  -  V( II) ) 

Remarks .  RATIO  will  be  used  in  the  interpolation  in  TAB2V  or  TAB3V 

Subroutine  OUTPUT 

Purpose,  prints  out  paged  output 

Usage,  call  OUTPUT  (IOPT) 

Description  of  parameters. 

IOPT  -  SWATH  sets  =  2  when  the  run  terminates.  If  IOPT  =  1,  the  pages 
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are  printed  out  even  if  not  full,  which  may  happen  if  run 
terminates.  Otherwise,  OUTPUT  does  not  print  until  a  full 
page  of  40  lines  is  reached. 

Remarks .  1)  Page  1  contains  a  time  history  of  some  ship  state  values 

and  the  position  in  fixed  space 

Page  2  contains  a  time  history  of  the  control  surface  values 
Page  3  contains  a  detailed  time  history  of  propeller  forces 
and  moments 

Page  4  contains  a  detailed  time  history  of  rudder  forces  and 
moments 

Page  5  contains  a  detailed  time  history  of  fin  forces  and 
moments 

2)  IOPTN  (2)  controls  pages  3-5.  If  1,  pages  are  printed.  If  0, 
they  are  not 

3)  Output  is  in  MKS  system  with  angles  in  degrees 

4)  Time  history  values  are  loaded  into  page  arrays  and  printed 
out  one  page  (40  lines)  at  a  time. 

5)  See  sample  run  for  page  listings 
Subroutine  PROFOR 

Purpose,  to  calculate  propeller  thrust,  pitch,  and  yaw  moments. 

Usage.  CALL  PROFOR 
Subroutines  Called.  TAB2V 
Remarks .  1)  Thrust  =  Kypn^ 

where  Kj.  =  thrust  coefficient  obtained  through  table 
look  up 

n  =  revolutions/sec 
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0  =  propeller  diameter 


2)  Kj  is  obtained  through  interpolation  using  TAB2V  and 
propeller  K^.  vs.  J  curves  as  data 

3)  Propeller  Y  and  2  forces  have  been  neglected. 

4)  Propellers  are  operated  independently  so  that  differential 
thrut  can  be  applied  for  low  speed  turns. 

Subroutine  RUDFOR 

Purpose,  to  calculate  the  rudder  forces  and  moments 

Usage.  Call  RUDFOR 

Subroutines  called.  TAB2V,  TAB3V 

Remarks .  1)  RUDFOR  first  calculates  the  aft  rudder  contribution,  and 

then  the  fixed  forward  rudder  contribution 

2)  The  aft  rudder  data  is  read  in  as  four  linear  arrays 
AXRUD,  AYRUD,  AKRUD,  ANRUD.  These  arrays  contain  rudder 
force  data  based. on  rudder  deflection  angle  and  speed. 

3)  TAB2V  performs  a  two-way  interpolation  on  the  rudder  data 
to  obtain  the  aft  rudder  X,  Y,  K,  and  N  forces.  The 

M  moment  is  calculated  as  the  X  force  times  the  moment 
arm. 

4)  For  the  fixed  forward  rudder,  an  array  BYRUD  has  been 
constructed  that  contains  the  rudder  Y  force  based  on 
speed,  drift  angle,  and  yaw  rate.  TAB3V  performs  a  3-way 
interpolation  to  obtain  the  rudder  Y  force.  The  X  force 
is  approximated  as  Y/7  and  K,  M,  N  moments  are  calculated 
using  the  X  and  Y  forces  times  the  appropriate  moment  arms 
This  portion  of  modeling  is  still  under  development. 


5)  Different  aft  rudder  sizes  are  handled  by  an  ARATIO  term 
in  the  aft  rudder  equations.  ARATIO  is  a  ratio  of  actual 


rudder  size  to  a  standard  size  rudder  used  in  the  model  tests. 
Subroutine  TAB2V 

Purpose,  performs  a  two-dimensional  interpolation  on  a  given  data  base. 
Usage.  CALL  TAB2V  (VI,  V2,  OB,  Nl,  N2,  ARGI,  ARG2,  ID,  KEY,  II,  JJ ,) 

VAL,  IER) 


Subroutines  Called.  NGHBR2 
Description  of  Parameters. 

VI,  V 2  -  variable  arrays  starting  with  smallest  and  increasing 

monotonically.  Values  must  be  unique 
OB  -  observed  function  values  stored  columnwise  according 

to  VI,  V2 


Nl,  N2 
ARGI,  ARG2 
ID 
KEY 

II 

JJ 

VAL 

IER 


number  of  Vi's  and  V2's 
arguments  1  and  2 
unit  for  printout 

if  =  0,  do  not  perform  search,  use  previous  V 
and  ratios 

closest  point  to  ARGI  S.T.  V ( I I ) . L  E.  ARGI 
closest  point  to  ARG2  S.T.  V(JJ).L  E.  ARG2 
interpolated  value 

0  -  in  range,  otherwise  out  of  range 


's 


determined 
by  NGHBR1 


Subroutine  TAB3V 

Purpose,  performs  a  3-dimensional  interpolation  on  a  given  data  base. 
Usage.  CALL  TAB3V  (VI,  V2,  V3,  OB,  Nl,  N2,  N3,  ARG2,  ARG2,  ARG3,  ID,  KEY, 


II,  JJ,  KK,  VAL,  IER) 
Subroutines  Called.  NGHBR2 


Description  of  parameters 


VI,  V 2,  V3 


OB 

Nl,  N2,  N3 
ARG1,  ARG2, 
ID 
KEY 

II 

JJ 

KK 

VAL 

IER 


-  Variable  data  arrays  starting  with  smallest 
and  increasing  monotorically.  Values  must  be 
unique. 

-  Observed  function  values  stored  columnwise 
according  to  VI,  V 2,  V 3 

-  number  of  Vi's,  V2  s,  V3's 
ARG3  -  arguments  1,  2,  and  3 

-  unit  for  printout 

-  if  =  0,  do  not  perform  search,  use  previous  v's 
and  ratios 

-  closest  point  in  VI  to  ARG1  S.T.  V(II).LE.  ARG1 

-  closest  point  in  V2  to  ARG2  S.T.  V(JJ).LE.  ARG2 

-  closest  point  in  V3  to  ARG3  S.T.  V ( KK) . LE .  ARG3 

-  Interpolated  value 

-  0,  in  range,  otherwise  out  of  range 


Subroutine  TOTVEL 

Purpose,  to  calculate  Inflow  velocity  at  a  given  point  based  on  the  radius 
vector  and  the  ship  velocities. 

Usage.  CALL  TOTVEL  (UU,  VV,  WW,  PP,  QQ,  RR,  XDV,  YDV,  ZDV,  UTOT,  VTOT,  WTOT) 
Description  of  Parameters. 


UU,  VV,  WW  -  ship  U,  V,  W  velocities 

PP,  QQ,  RR  -  ship  P,  Q,  R  angular  velocities  (radius/sec) 

SDV,  YDV,  ZDV  -  X,  Y,  Z  coordinates  of  point  in  question 
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Remarks. 


UTOT,  VTOT,  WTOT  -  U,  V,  W  inflow  velocities  at  that  point 

1)  All  velocities  are  in  the  body  system 

2)  UTOT  =  UU  -  RR-YDV  +  QQ-ZDV 

VTOT  =  VV  -  PP-ZDV  +  RR-XDV 

WTOT  =  WW  -  QQ-XDV  +  PP-YDV 


3 

i 


1 


J 
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RECOMMENDATIONS 
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The  current  digital  simulation  can  be  used  to  develop  optimum 
maneuvering  strategies  for  various  missions  to  be  performed  by  SWATH.  The 
results  of  such  studies  if  they  are  conducted  may  provide  useful  information 
to  the  Navy  allowing  a  complete  evaluation  of  SWATH  potential  capabilities. 

Another  major  area  of  mathematical  modeling  that  should  be  pursued 
is  the  modeling  of  sea  way  effects.  Due  to  the  small  waterplane  area,  it 
is  expected  that  seaway  effects  would  have  a  relatively  mild  influence  on 
craft  dynamics  as  compared  with  conventional  craft.  Taking  these  effects 
into  account  would  lead  to  a  more  meaningful  comparison  of  SWATH  performance 
with  conventional  craft  performance. 

There  are  many  minor  reffnements  that  can  be  made  regarding  the 
mathematical  model.  The  mathematical  model  and  digital  simulation  program 
are  designed  to  accept  refined  data  bases  and  improved  theoretical  estimates 
when  they  become  available. 

It  is  recommended  that  the  basic  data  base  approach  employed  in  this 
study  be  extended  to  other  high  performance  craft  that  are  currently  employ¬ 
ing  theoretical  models.  These  theoretical  models  are  usually  very  complex 
and  require  extensive  computational  efforts.  However,  one  could  limit  these 
computational  efforts  by  using  the  theoretical  expressions  offline  to  gener¬ 
ate  a  data  base,  a  one  time  effort,  that  could  then  be  used  in  the 
approach  presented  in  this  paper  to  carry  out  aany  different  simulated 
maneuvers. 


CONCLUSIONS 


j 


This  study  represents  a  fairly  extensive  dynamic  simulation 
model  development  based  primarily  on  experimental  model  data.  The 


82 


final  product  has  proved  to  be  a  useful  evaluation  tool  in  studying 
SWATH  performance  during  a  variety  of  maneuvers  for  several  different 
configurations.  The  success  of  the  approach  used  is  based  on  a  well 
planned  program  of  model  experiments  that  provided  a  sufficient  data 
base  for  mathematical  model  development. 

The  results  of  the  predictions  indicate  that  significant 
improvements  in  performance  can  be  achieved  when  the  spade  rudder 
configuration  are  considered  or  when  a  deeper  draft  is  used  during  turn 
maneuvers.  The  results  also  indicate  that  the  SWATH  has  very  good 
turning  performance  characteristics  at  low  speeds  when  differential 
thrust  is  used. 
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NOTATION 


S  State  vector 

A  Generalized  mass/inertia  matrix 


F  Time  rate  of  change  of  state  vector 

M  Ship  mass 


X. 
u 

Y. 

V 

Coefficient 

used 

in 

representing 

X 

as 

a 

function 

of 

Coefficient 

used 

in 

representing 

Y 

as 

a 

function 

of 

Z. 

w 

Coefficient 

used 

in 

representing 

Z 

as 

a 

function 

of 

K. 

P 

Coefficient 

used 

in 

representing 

K 

as 

a 

function 

of 

M. 

q 

Coefficient 

used 

in 

representing 

M 

as 

a 

function 

of 

N. 

r 

Coefficient 

used 

in 

representing 

N 

as 

a 

function 

of 

XG 

Location  of 

CG  along  body  x  axis 

yg 

Location  of 

CG  along  body  y  axis 

ZG 

Location  of 

CG  along  body  z  axis 

Y. 

P 

Coefficient 

used 

in 

representing 

Y 

as 

a 

function 

of 

Y. 

r 

Coefficient 

used 

in 

representing 

Y 

as 

a 

function 

of 

z< 

Coefficient 

used 

in 

representing 

Z 

as 

a 

function 

of 

K- 

V 

Coefficient 

used 

in 

representing 

K 

as 

a 

function 

of 

K. 

r 

Coefficient 

used 

in 

representing 

K 

as 

a 

function 

of 

M. 

w 

Coefficient 

used 

in 

representing 

M 

as 

a 

function 

of 

N. 

V 

Coefficient 

used 

in 

representing 

N 

as 

a 

function 

of 

N. 

p 

Coefficient 

used 

in 

representing 

N 

as 

a 

function 

of 

I 

XX 

Moment  of  inertia 

about  x  axis. 

‘yy 

Moment  of  inertia 

about  y  axis 

>zz 

Moment  of  inertia 

about  z  axis 

u 

v 

w 

p 

q 

r 


P 

r 

q 

v 

r 

w 

v 

P 
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Product  of  inertia  about  xy  axis 
Product  of  inertia  about  xz  axis 


Product  of  inertia  about  yz  axis 
Angular  velocity  about  x  axis 
Angular  velocity  about  y  axis 
Angular  velocity  about  z  axis 
Velocity  along  body  x  axis 
Velocity  along  body  y  axis 
Velocity  along  body  z  axis 
Roll  angle 
Pitch  angle 
Yaw  angle 

Displacement  along  fixed  x  axis 
Displacement  along  fixed  y  axis 
Displacement  along  fixed  z  axis 

Transformation  matrix  from  body  to  fixed  coordinate  system 

Effective  horsepower 
Total  force  along  body  x  axis 
Total  force  along  body  y  axis 
Total  force  along  body  z  axis 
Total  moment  about  body  x  axis 
Total  moment  about  body  y  axis 
Total  moment  about  body  z  axis 
Density  of  water 
Ship  length 


F 

AR 

Vf 

CL 

a 

A 

X, 


f 


D 

J 

V 

Z 

z 

M 

M 


a 

9 

z 

9 

z 


Fin  force  or  moment 

Fin  aspect  ratio 

Inflow  velocity  at  fin 

Fin  lift  coefficient 

Fin  angle  of  attach 

Fin  area 

Fin  x  position 

Fin  y  position 

Fin  z  position 

Thrust  coefficient 

Propeller  revolutions  per  second 

Propeller  diameter 

Advance  coefficient 

Propeller  speed  of  advance 

Coefficient  tr  represent  buoyancy  Z  as  a  function  of  9 

Coefficient  to  represent  buoyancy  Z  as  a  function  of  z 

Coefficient  to  represent  buoyancy  Z  moment  as  a  function  of  9 

Coefficient  to  represent  buoyancy  Z  moment  as  a  function  of  z 
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Fin  force  or  moment 

Fin  aspect  ratio 

Inflow  velocity  at  fin 

Fin  lift  coefficient 

Fin  angle  of  attach 

Fin  area 

Fin  x  position 

Fin  y  position 

Fin  z  position 

Thrust  coefficient 

Propeller  revolutions  per  second 

Propeller  diameter 

Advance  coefficient 

Propeller  speed  of  advance 

Coefficient  to  represent  buoyancy  Z  as  a  function  of  0 

Coefficient  to  represent  buoyancy  Z  as  a  function  of  z 

Coefficient  to  represent  buoyancy  Z  moment  as  a  function  of  0 

Coefficient  to  represent  buoyancy  Z  moment  as  a  function  of  z 


APPENDIX  A 


DATA  BASES  USED  IN  THE  SIMULATION  MODEL 
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ooon  o  o  o  on  nnr>  noon  oo 


E°S  f  1  -1  > 


-  p  rL  A  T I VE  ERROR  CRTTEPIA  pC®  Y'l-IE) 


COMMON  'xccntl/  NCONT.I 
t  CV2< 10) 


COMMON  /  CU  T  X 


ONT.IC1  (10  )  ,  IC2  (10  )  ,  I  A^Tl  1  r.  >  ,  JACTfio  >  ,CVlf  10) 
2« 10) ,CV3( 10),CV4( IQ) ,CPTNl ( 16) ,CR IN2( 16), 

IN  1 (1 ?) »C0TN2(12)»C0If3<12) ,°ELAB ( 2) »CHFC*<2> 
TE » 12 > .KACT( 12  ) 


COIN  1  (l?)  ,C0!N2(  12)  ,001*  3  (12)  ,or.U 
PATE ( 12 >  .KACT( 12 ) 

IBU3  (5)  ,TPP  I  NT  (5  )  ,  JR  A  ,*|R  t**  ,N®  R  TNT 


COMMON 

niMENSI 

F 

A 

define 

U  - 


R  - 
OMJ  s 
THE  T  A  = 
®SI  = 


TIME=T 


/HEADER  /  TITLr(10) 

ON  Y  ( 13  ) »A(6,6)»E(13) 

-  GENERALIZED  FORCE  VPCT0R 

-  GENERALIZED  M A SS/I NE« T I A  MATRIX 
LOCAL  VARIABLES 

Y(?  1  09 

Y(3  )  USli 

Y  (4  ) 

Y(5  ) 

Y  (6  ) 

Y  (7  ) 

Y(B  ) 

Y  (9  ) 

Y  <  1  0  ) 

Y  (11) 

Y  (  1 2  ) 

Y  ( 13  ) 


DRAFT 


COMPUTE  S°EED  AND  OTHER  REQUIRED  VALUFS 

U2  =  U*U  ♦  V*V  ♦  W*W 
SPEEC=  SGRTfL?) 

D RIF  AN  =  -ASIN(V/SPEEO) 

COS°HI  ^  COS ( PH  I ) 

COSTHE  =  COS (THETA ) 

COSPSI  =  CCS(PSI) 

SINPHI  =  Sit!  (PHI) 

SINTHE  =  SIN (THETA ) 

S I N  °  S I  =  SIN(PSI) 

COMPUTE  FIXED  TO  BODY  TRANSFORMATION  MATOIX 


tmtrxc  i, i) 

TMTR  X (  1,2) 
TMTRXC1,?) 


COSTHE  *C0  SPS I 
COSTHr*SIMPSI 
■SINTHE 


TMTRx (2,1)  =-C0SPHl *SIN®SI+  S INTHE*SINDHI * CO SPS  I 
TMTRX(2,2)  =  COSPHI *C0CPSI  ♦  SINTHE *SI NPHI*S INPSI 
TMT  R  Y ( 2 , 3 )  =  COSTHE*SI WPHT  ' 


TMTRX (3,1) 
TMTRX(3*2) 
TMT  R  X  (  3 ♦  3) 


SINPHI *SI NPSI 
■SINPHI *C0SPST 
COSTHE  *C0 SPHI 


SIMTHE*COScHT*COSPSI 
SI  N TH E  * C 0 S P H I  * S I N 9 S I 


COMPUTE  CONTROLLER  ENVIRONMENTAL  cACTORS  SUCH  AS  WIND  OUSTS, 
COMPUTE  PILOT  ACTIONS,  major  LINK  IS  VIA  COvMON 


on  o  o  o  o  o  o  »-•  o  r)  o  n  on  o  ">  n  n  o  o  o  n  o  non  o  r>  n  o  n  o 


CALL  COVTPL(T*Y> 

COWDUTE  TOTAL  FIN  FCRCE^  AND  «n^‘'MTS  . 

CALL  FIN^OR 

rn,MP-jTE  TOTAL  RUDDER  FORCES  AND  MOMENTS. 

CALL  R  UO  cO  R 

f 

COMPUTE  RARE  HULL  FORCE1;  AND  FOMENTS  ; 

CALL  HUL^OR  j- 

COMPUTE  AERODYNAMIC  FORCES  AND  VO^^NTS 

» 

call  AIRc0R 

COMPUTE  2 U C Y A  N C Y  FORCES  AMD  M0MrNTS  ■ 

CALL  3  D vFO  R  j 

COMPUTE  cR0PFLLE°  FORCES  AND  MOMENTS  ■ 

CALL  Ponc0R 

SQM  P  L'TF  ACCELERATION  AN”  ADDITIONAL  -ATE  DEPENDENT  FORCES  AND  MOME 
CALL  ACOFCQ  J 

SUM  TOTAL  FORCES  ANC  MOMENTS  \ 

CO  1  1=  1*  G 

r<I>  =pFmcI)+«MJDFM(I)+HULLFMCI)*AI»pm<  I)^PUOYPM(I) 

4  *DR  OP  PM  (  I )  ♦  AOOFMCT) 

CONTINUE  l 


A  OD  ON  TrR  MO  DUE  TO  CENTER  CF  GRAVITY  NOT  AT  ORIGIN 


F(1):FU) 

F <2  >=FC2> 

F<T  )=F  (  T  ) 
F(R)=FfA) 
r  f  5  ) =F  Cc  ) 

F  <S  >  =F  ) 

ADO  ON  ADDITIONAL 


♦  MASS*  (XC*(Q*Q-»'R*D>-Yr*p*(3-ZC*P*R  > 

♦  MASS*<YC*(R*P+P*C5)-ZC*a*R-XC*Q*P) 
♦«ASS*(ZC*rP*P+Q*0)-XC*R*P-YC*R*G) 
♦MASS*C?C*( U*R-U*P)+Yr*( U*G-V*P) ) 

♦  wASS*fXC*(  V*P— IJ*G)*-ZC*( V*R-W*Q) ) 
+MftSS*<YC*< W*C-V*P) +XC*( v*p-u*° ) > 

TE°  MS  ARISING  FROM  SPACE/EOOY  TRANSFORMATION 


F(1  )  = 
-C2»  = 
F(3)  = 


C(1>-KV*P 
pf2 >+(W*P 

c  m+cuo 


-  y*Q)  ♦■'ass 

-  U*R) ♦MASS 

-  V*P) *MASS 


p<A>  =  ( T  YY  — I ZZ) *0  *9  ♦ 
R<5>  =  C TZ2-IXX) *R *°  * 
Ff«.)  =  <TXX-!YY)*°*Q  ♦ 


IXZ*P*a+  F<*>*  fO*f-R^R  )  M  YZ  -  ixy*p*r 
IXY*Q*R*  F(C)+  (R*P_P*P)*IXZ  -  IY  Z*G*P 
IYZ *R *P ♦  C(G)+  (P*P-Q*Q)*TXY  -  TX2*D*Q 


LOAD  INERTIA/*' ASS  “ATOIX 


'->H&  ’***•:&& 


•sm 

a' 


^  J*  MW 


9? 


MP  ""I'JJ 1 

-v>r'  -cV" 

•jf,  _  ,  •- 


“V" 


UUM.j.,  ii  n» 


rr 


V'r-y>^ 
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DismFT 

A(I,^)=  3.0 
ECN  TINUr 

A  ( 1  *  1  ) =M ASS-XUDDT  * 

a n  * c) =”srs*?c 

A  (1  ,c)=-i-a<;s*yc 
A<2  *2>=M*SS-YVDOT 

4(2*4)  =—P'ASS*ZC  —  YPDC'T 

S (2  *  4) r-VRPOT  +  VA SS  *XC 
4(3*3)  =J*aSS-?UDOT 
4(3  * 4 )  =Y ASS*YC 

4(3*5)=-wASS*XC-7QDCT 

4(4  *2)  =  ->‘4SS*Zr-KVJ0T 

4(4*3)=Y4SS*YC 

4 ,4) =1 yx-KPDOT 

4(4,5>=-IXv 

S(4  •*)  =-ry7-*°riT 

4<3*l)=‘-,ASS*?r 

4 (5  *3)=-vASS*XC-vWD0T 

;  (3.4)=-TXY 

5  (5  *5)=lYY-M0f)0T 

s  (5  *6)  =-TYZ 

4(6  * 1) =-*ASS*YC 

A(6*2)=  +  wASS*XC-4’VDCT 
i(g, 4)=-TX7-NPnoT 
4 (6 ,3) =-IYZ 
A ( &  * 6 ) =IZZ-NPDCT 
IF(ID9F.NE.e>GC  to  10C 
S0LV3  FOR  (DOT (  VECTOR 

CALL  LINQ<A,F,e,IE9R> 

IF(  IERP*’ E.O  )  GO  TO  1900 

DHI  DOT .  THETADOT*  PSIDOT*  IJF  *  Vr*'4E*  SPEED 

F(7)  =  (D*STNTHE*SI'<PHI  *R  *  S  IMTH  E*CO  SP  *  T  ♦  P*COSTHE>  /COSTHE 

F(3)  =  C*C0SPHI  -  P  *SI  4.  PHI 

F(9)  =  ( E*COSPH I  ♦  G*SIV°HI ) /COSTAE 


mg 

P 

||| 

12 

r#$ 

DO  30  1=1*3 
SS  -  0.0 
00  2S  0  =  1*3 

SS  =  SS  +  TMTPX(J*I)*Y(J) 

CONTINUE 
FC9+I)  =  SS 
CONTINUE 

=  (13)  =  S0FT(p(10)**2*F  (ll)**2-*-r(12)  **2) 

SO  TO  9a9 
COMTINLE 

RES  REE  S  °F  FPEEDOF  ARE  LESS  THAN  6 


s  -  *- 


U  O  U  U  U  C  —I  —  U  L.'  K  C) 


3G 


PITTN  AN  OLE  pIxFO.  THDOTz  0 .0  T^PLIrc  0  =G*COSPu  I-»  *  S  INPuj  •  HFNC® 
Cr>$onT*r.rnT-SIfs:PHI  *  =  ^OT-PHTOOT*(P,*.<;r‘lPHI  +  R*COSc>uI)»  °H  I''!}  T  TS 
USED  Fr  CONVENIENCE  •  EQUIVALENT  NON  CO^K  FX°BESSION  COULD  °E 
USE C  FCP  DWIPO  T. 


A(5«l)=r'.C 
S  (5 , 0> =3.0 
ft  U  ,<=)  =0  .0 
A(3  .5 ) =C  .0 
f  f4  *E)=n  .0 
a(S.F) =" .0 
a  (5*3>=C.O 
Af5,4»=0.3 

.  c  1 


DRAFT 


A  (5  ,A>  =-SIN»H 
comp  ut  r  FH  IOC  T 

e  <7  )  =  f Q *S  IN  THE*S IMP H I  +  P •*  S  TNTH F*CQ  CPH  I 
-  (5  )  =  c-f  7)*<Q*SP!PHI^P*C0SPWI> 

IFf TC0F.LT.5J  TO  ''■OO 
CALL  LI*'0<  A.F,6,IERP  > 

IF(  I  FRC.NF.01  PC  TO  l=r5  = 

C(S1  =  3  . 

"(?)  =  (R*COSOHI  +  3  *SIVPHU  /CCSTH^ 

CO  1 30  I  =  1,^ 
so  =  o.a 
DO  130  J  =  1,3 

ss  =  ss  ♦  T«Tpyr j,r )*y(-j> 

CO  NT  IN  IF 

c(a-*-n  =  ss 

CON  T INLr 

Ffl3)  r  SORTCF(lC)**2  +  F(ll)**?-*-Ffl2>**2) 


c*  COS THF ) /COSTHF 


~.n 


’ n  cqo 


CONTINUE 


IF( IOOF.LT .4 )  00  TO  40G 


ft  DOP  NO  PrTrH,  NO  UFA  VF  *  I','»LIrS  TH^TA  DO  T  =  3  ,  ,  WE  =0  .  . 
H=-U  *S  IvTUF+V*-COSTHr*SIf'PHI  ♦  W*COSTmE*CCSPPI  TmPLTFc 
UDOT  *S  INTHE  -  VD  CT*  COST  HE  *S  I  NP*-*I  -yDO  T*  COSTHE  *CO  SpH  I  = 
PHIDOT*V*COSTHE*  COSpHI-PHIFQT*U*COSTPF*SINPH  T 


'JFDOTrQ.O. 


ft  (3  »  1  )  =  FIMHE 
A <3  *2) =-COSTHF*SINPHI 
ft  (3 ,3) =-COCTHE  *CCSPH I 
A (3 ,5 ) =  0.0 
ft  (3, 4)  =0.0 
ft  f  3  ,  F  >  =3  .0 
? (4 ,3) =0 .0 
a  '5  =n  ,o 


CC3) =  Ff  7) *  CCS THF*  CV *COSDHI— W*SINDUI  ) 


CALL  LINPf  «  ,F,<,,IERR> 

IFf  IFP°.ft.F.0>  DO  TO  10°- 
F  («  )  =  0.0 
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r>  t-* 


=  O  )  =  (  "  *C  E I  D  *SI  »'PH  I }  'Cf'^  Thf 

10  3  30  T=1*T 
SS  =  0.0 
DO  3?0  J  =  l»* 

SS  =  CS  •►  T"*T^x(j,T>*y(j) 

3  •'9  ~r.MTI'i;jr 

r  (  ~  ♦  I  >  =  0  C 
3  3  0  COM  T  TN LE 

-(12)  -  ".o 

F  (1  3  )  =  SGRT(F(  10)  +*?*F  (  1  1  )**?) 

0,0  TO  09° 

400  rCNTIMLT 

W  D I  T  E  (  6  »  6  0  1 )  IL'OF 

6  01  P0R  M  M  (2  4H  VCEEL  MOT  AVAIL5qLE  r0!7*l4tlrtJ  OF  FR  E  E  DO  U  . // ) 

C 

cTiP  6 
C  C 

=  09  C ON  T  INUE  I; 

I  c<°  NT  (!•:>>  00  TC  10  "*C 
p  ET  U  c  N 

c  ***  I  McpMFO  1ATF  OEPUD  "UTSUT  ***  *** 

r 

1000  '•iRITF(TD13*6r,3)  T,U*V,W,PHI.THFT«.D9TFAN,  (  CF*.(  "  )  *°  IJDF*  (  *  >  ,HMLL~ 
1 (K) , £ I  RE  MR ) ,BUOYFM(K)t * TCF* ( K ) * PROP E« ( K ) , K =1 *0) 

PNT  t  10  >=. FALSE. 

RETURN 

0°p  CONTINLE 

\'0  INVERSE  F CR  MASS/TNER-IA  * A TO  IX. 

’JRITE<  IPEtOOO) 

'JRITr(  I®  6,  £C?)  (FCI),I  =  l,6),((flfTt.J)  .0=?  ,6  >,1  =  1,6) 

W»ITE(  T?F,?0?)  rvx,IYY,I7Z»tXY,iy2,r  YZ,n'COT.YVOrT.YDDnT,Y?.r  "T* 

1  ZWDOT.ZOOCTtKVnOT,KPDOT,KTDOT,3yDOT,MODOTtMVnCT, 

7  MPDOT,MRQOT 

6 "  0  ^ORYSTtS'H  I?,VE9?E  OF  V  4  OS  / 1  NEDT  I  A  VMTETV  DOES  NOT  EXIST.  ETOP/) 

602  FOR^ATtEP  EALX,1CE12.4) 

6  03  FOR^?T(T?*4HTIME*T10*lHL'.Tl6»luV«T22»lHy.T27»3HPHI«T32» 

S9HTHE Ta.T3  9tAH?ETA,TElt7MFTN,T5r,*6HPUDDEr  ,T71  »4HHULL, 

PT82  »4HAEPC*T=9,8HBUOYANCY.'r105«4HRATE.Tl  1  3  ♦  1  OHPR OPULE T0 U  ./. 

C1H  ,  F5.1  »c>Eb.2*T45*lt-Y»  7E11 .4,/,T4EtlHY.7El]  .4./, 

CT45»1H2.7E11  .4. / tT4E * 1H * *7E 1 T .4. / ,T45 ♦ 1 p* *7E I  1.4*  /  *  T  4  r  ,  i lJ  N  *  7  E 1 1  . 
C*/) 

STOP  5 
END 


PUMCTIC!"  DPAACG(V)  A#lft/lf  J 

DTRCOR^r  fTEPPOLATI  "N  ?r0  FOR  SWATH 

n  t  \*  r  ».i<5  ^  f  \  RPrEC  f  25  )  «  FHP  (  Rc  ) 

CAT  A  SPED 

.  1 «?.  ,20.*21  .  »22  .*2T  .«24.  ,25,  f2F  .  »27.«2R.  t 

-iTA  F{-r/O..ir>C..2OC.,AnC.ta00..1200..1Fr0.»?-5R9.,S^0..^n0., 

J  "  5  n2  2  .  tSc  o  C .  ,  6  0  0  a  .  ♦  €  7  * o . *7« 0  0. • Q50  ? . . 1 2 0 0  C . ♦ 15  2  0  0  . » IP h 0  0  .  * 

1  21900..24  7CO.*27&CO.*3COOG.«3  2f?OQ.*3  7POO./ 

Vw  s  =v 

'/=V>*S*1.F44 

IF( V.LT.C.)Or  TO  99 

00  10  1  =  2 1 25 

ICC  V  .LE. SPEED ( I) )G0  TO  ! c 
to  rOMTINLr 

15  't  AT  I  C  =  fV-SFEED(I-l >>/<SdEED< I) -SpEEDt I-l> > 

rMpv  =  ( CHP ( I)-EHP( T-l)  >*R AT  IOfFHPf I -1 ) 
dp.aa^g  =  (FHPV*145Q.  1  >/v 
V  =V  **  C 

p  £T  U  R  *! 

Q9  USITr(6»20>V 

20  =  QR  '*  ATI*  ?  F c  E  0  EXCEED  S  RANGE  *»F10.41 

RETURN 
END 


o  r>  n  ri  o  n  n  n  n 


rUPP  CUTTNE  ELLFP  (T,  °t,y  »KpTATn  ) 
rUL  EP  I^TETP ATION  ROUTT  NT. 

DIMENSION.  Y(NSTATF),C(^''> 
rnMoijTc-  chs  cc  differential  F0UATI0»'c. 

CALL  TAUy<T,vtF> 

TNT  FOR  AT  r  FROM  T  TO  T  ♦  DT 

00  1  I=1,NFTATE 

Yf I)  =  VCI)  ♦  F(I)*nT 
1  CONTI  M  L  F 

C 

C  UPDATE  TTWF 

C 

T  =  T  CT 

RETURN 

END 
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DRAFT 


SUP  o  OU  7  T  \p  FINFCR 

cco^cftM  TO  CCMPUTF  57 1 'C  cCOCFp;  ftNO  MOUNTS 
CQM  M  ?N  /CONST/  <»ft0i;vi*»H'‘A,RF0ft?,CH0  V  ,  °  H^W ?  ,  r  TTT1 
~GMMrn  /FT4TE/  U.V,W«P,0*c.P-JI,THFT\*^F  T»XF*vF*Zr»FPrEC*<J»0°.  IFftN  , 
lFOS^HI  ^CCSTHE^CO^PSI  *SI  VPHI  tSI*jTHE,S  IN'P0!  *T‘  TRVf  .TIMt 
COMMON  / c I  NV  i°  /  XFIMC*)«YPI**f  *1  *?CI*  t *  ?  •  A*  A®  t *1/ 4  )  * CLclS  <4  I  •  0?F!  *1 
\  <  4)  ♦  FT  *>  A  'G  (  4  )  ♦XF2RF  \i  (4)  ♦  YP09FN  (  4  )  ,JF  "3FNf  4  >  ♦  A  LIP  T(  4  )  ♦r'»  ft  0  (  4  )  .PPM 

!  (F)  ,cif*.  \n<?> 

no  i  »i=i  ,4 

CALL  TCT  VFL  (  U  •  V  ♦  W*P  *  Q,P  »  XF  J  N  ( *M  ,  yp  I  *J  (  V  )  ,  ?F  I {  N  »  .  'JTOT  a  L  ♦  V  TO  T  AL  » 
WTOTAL  5 

GAM  =  ATA»’  fWTOTAL/UTOT  AL) 

ALFA=  GAM+FINA  NSU'» 

TRF  0  R  =  GAM 

ALIPT(  N)  =RF0W2*  ABAFT  M(N  )  •  (U^OT  ♦  ’.♦,*t''T!>l»*2.>*CLp:N{N>*ALFA 
DS«G<M)=aBC(ALI=T(N))/7. 

XFOR  FN  <\  >  =  AL  IFT  (  V' )  *0  IF'  (  TP  F"  ®  )  — CP  A  0  (  *J  '  *  r  (  Tc  c  T  3  I 
7POPF'4<N>  =  -CALIPT(N)*COr(Toc''’p)*nPAr,  tm*-;*  (T 9® o 3  )  > 

1  CON  TIN  LE 
00  2  N=! ♦  F 

2  PP”  f  *>)  =  0  • 

0"  3  N=l*4 

PF5*  (  1 )  =F?“Cl  )-*-XPOPP\  (Ml 

3  FFM  (0)  :Fcw  <  3  WFORFN  (N) 

on  4  V  =1  ,4 

CFM  f  4)  =  pp“{4  )>ZF  CRFN  <N)  *  YP  I M  f  *J )  -  Yr0°  F  *J  {  *  )  »Zr  I '!  (  N  ) 

FFM  f  P  )  =F  =  W(F»-*-XFC®FA4  <N)  *7FIN  (*J  )-7P0R  CN  (F  )  *YrI'J(N  ) 

4  FFM  (  F  )  =FF«  (F  >+YFCRF'l  f  N)  *YFIN  <N  )-YFOR  C*J  <N>  *Y®I  N('J  ) 

PETUPM 

END 


o  r> 


SIJRR  GUTTf.T  HlLFOp 
°EA  L  LENGTH 

REAL  MAC.  S*IXX*TXY*IXZ*rYY*TY7*IZZ 


DRAFT 


COMM  ON  /  T  T  A  T  E  /  U»V *U«P «O.R *RH T*THFTA*PS  I  *XE «YE, ZE*$PEEQ *E* 

1  rRIFAN,CCSc*HI,CPE;THE:,COSPSrf  CrNPHI,^TNTHF,SINP?  T  , 


T«TPXf3»T),TI*r 

COMMON  /rO'-C  7/  SADIA?'*  rPOC,  BHriA2»PwOVi,PP3w?,FTTOM 
common  /rpwnvF/yc*YC*?c 

COMMON  /VA9"Lt/tF  A)*  DUM*»V(44) 


COMM  CM /H  VCP0/LGC4)  *L7*KY(4,321  3)  ,c{4  ,  \c  )  .x^f  a,->q  )*  TPCINT  (A)  , 
1U.MA  X ,HYPSFM<  ft) ,ICNT 

POYMON /CP ART /SC FACT  *LFN PTH ,M » SS* I  XX « I  X Y  * TX Z . T YY t IYZ *  I II 
COMMON  /HCLVAP/HMLLFMf 6) * ^UL ° I K 
"O  45  K!J  =  2*g 
45  MJU.FM  <K»j>=o«3og. 

H(1  )=ppr*9  40  TAM 
H(2  )  =OPIFAN  *RACT  AN 
h  f3  >  — R  ♦LFNGTP/L1 


,°<l4 

CALL  HYHIMT 
r.SS=SQPT^J*U+V*V) 

CCF  CEr-O^AASG (SSS) 

HULLFM  ci  )  =  FOPCF*a  *«•  (.  0  ?*  CHIJLP  l*  *  7E>  ) ) 

'-ULLFM  <2>=HYC*FN,f2> 

u  UL  L  F  m  (7  I  =  C  #5**1?5*U)  *  R  F0W2  *LEN5  TH*L  ENG  TP*  .  0  C  1*U*U 
HUL  LFM  f 4  ^  =  FYCRF  M  f 3 ) 

FULLFM  (5>=HULLF”U)*f  R-^F-C  (  UULP  IX+ZEI/S*  )  ) 

-ULLFM  ff  )=R YCRFM (4) 


T5RM$  C LE  TO  CENTER  OF  GRAVITY  NOT  AT  ORIGIN 

HULL  FM  f4)=HULLFMI4)  +  MASS*°.,«1  *  /  YC *T.M TR  X  <  3  ,  .3  >  - ZC*  TMTR  X  ( 2  ♦  7  >  ) 
M  ULLFM  (-  1=H1JLLFM  (5  )  -v  AS  S*o  .R  1  *  (X  C*TMTR  X  /  3 , 3  ) -ZC-TMTR  X  fl  ,7)  > 
H ULLFM  (c  )=HULLFM(6)*WASS*9.S1*(XC*TVTRX(2*3)-*YC*™TRX(1  ,3)  ) 
°ETURN 
TND 


r 
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r!-~~K  .  T- C!“£L  *?P  "«T  •  !  e  TL'"  TW  ') 

1  ’■  l ir  s-T  tc  & l v r *■ T  '•sy  wcl,,i:.<t!'  »*a  -*ta 

i  a  <t  -  (  r  >  a  > 

T~r  ■  f  6  1  i  •-'^/ii  =  y-r,  <  a. .  l  f-4  i .  r  T  ^  ?  r>  i 

.......  c-TrK  T c  Clr'~i‘"r  Viur  I' prifjiir-'T  V :  °  !  i  ‘  L'l':  (-<) 

•  :t  I’-  tmF  ,'i;nc  c  ',T  r^Lt-  (  'l  -  "  I  "  ~l  F  T  u  T  **  L  I'‘Tr^'/*L  1 

jr  *,rr'  -j  v  xi «  |_nrp  ,-tr  rnpt v  ▼  #  tc  »  r  t  PT,n  fry  PTLL. 

f  A  ,  — 

'  ■»  ?  '  ’ri.'i 
~  -H  (  I ) 

•  =T  ~  ‘  T  r. '  •'  T  ' 

I  r  <  T  r  ic 

5 ,T  (  x-.<T  >  ♦~T.2)  ’  .1  ~  ^ 

crjsri!  ijnyior  fc  r  t  •  a  ijry  iyp  j->  or*  T'iT» 

’  c  -  _  i 

*-y  » i 

'  c  y  : T  r  ?  '  =  ?' 

t- f  t  tj  >  :fo:vt(i  >o--i 

Tr( "  *f  ■  •  i.-  fT>>  t  '  :r 

^  '■  T  ->  1  - 

rr^cu  rov'<yaen  ^  -T-n  .vry  -,cn  ^^mt,  . 


7  1 


—v  _  1 

T  OA  T  ■  f  ;  T  ) -« 

:-*(  jn  ? o 

a  -T-  n  i 

rCMTI'-U  r 

.-TT"(  ■  .  1 1  "  '  (  T^H  ,'  =  ]♦'*) 

-  C  3  ••  •  T  '  *  -■ r'  I  V  T  FP  *.4TF) 

-  -  -  -  -  -  if  isi ,' =  •>  >■['.'  cr*rr ” Tr.*" v s  vi  st  c^Meu 

frc  I  ^  !  >  r.’LL 


~  a  7  -  -  -  i  ,  a 

■-  - 1  *~  ■'  I  -  -  :  I ) 


:?•  ‘"'UT-  ■?!  *T£-CF  T~  Pc  !  r  ;f  7. 
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>-<  ct 


DRAFT 


c  | 1  r>  r?  r  1 1  ' 


r  r  r  r  r  t  c  4  •  t 


ci  )'’i)TT‘  f  T"  r . ir  L  T  r  ~  D  !vTrr=',1t-'TI*‘- 

_  2.  ,  •  c,,,rf  -  p  ^  t  •!  -r  p  ■-*  r  ▼  '■  f  D  J  *  (_r  3  •  T""CTT 

■  *  “  <-  T  7  3  '  r  I*'  A  T  *  f '  L  £  . 


..  ..  vivf  f  i.Li  »LI.A  .!  L  :r  ~  •  -v  '  u  !<*•  ’■'->•  y  ~< 

cj  '  ;  ■ •  t  '  /.  '  ,  W  >  ,  -  vr~  v  ■  '(4'tir  T 

•  T  /  r-  ■  c  T  •  I  '  (I  ,  U  ,  ,  <  *0 


_  . PI-.'TI-'  T'  rr'-'"  JT-  ru'  ~  C?.  T  r  T 

•/»U  U!  13-1  "l  ?  *  <  T  v -1  'LL  T  *  ?  I  ^  - 1  >  '  5 
- (- »  t  c  *  LL*'  Th  ir  '•  r u  ” T '• r . 

T r i c ’ : ' + 1 

•;'!  T  ■  (  ,  7  : )  ^ 

'*-<rT  F”  CT  T  ”  V'-UT^  *T  cr- r  3"’  ~  T':~  3?I~  Df,TfJT'' 

:  p  ?  '  t  *  =  i  •  ? 

•']:  I  r  ?  T  I  "  f  *  '  *  T  7  -  t 

,  „  „ .  7 ' = i .  r 

j  2  T 1  =  1 »  ? 

■  x-  T  -  r  T  -'•»?  +  :  i  -  l 

? -  i:=i.° 

✓  a-  '  3  ~  T  : ;  T  (  t  )  -»  I  C  -  1  .  . .  ,  „ 

i/ 1 ,  i/o  ,  /  >,v/i,v.(i.i/i  ) ,  v  -  i"' ,  k  -  > .  '•  r-  (  T  . u' T  ) «  v  <  4  •  <  4  1 


Tf  ,  .  <•  ij  r.  <-  rq  j  •_  t  r  **4TT*B:?*  *  5  fM..I  7r !.  =  *  •  c  > 0  •  4  • 

1  *  -v.  T  e  7  7FL=+«Fl?«*t  *  T  •*.  *  «r*.r 13.4.  »  S‘  t  -  n  =  *  •  ‘  1  r  • ' 


!.  1  =  !  3  i  1  *  1  »  J  ?  •  ^  3  •  ^ 4  1 
i  =  i  *  1 
;  3  !7  =  i 

Tr t  Ll.Lr.r:  '?C  TC  1 r 
'  'I  "  1  >  -M  v  f  T  v  *  L  1  ) 

■  '  T  '  1  . 


.......  VALUES  en»  -F.*C  3r3Er>  .if'E  SF  T  rr  T  ST'?*' 

:(!'/.  •  n  =  o . 

“rvir*;'." 

Tr'E.’1'"  Mi •  li *  f  v-,»i  >  ,T  =  i  .4> 

"  AT  t  »  '  1  =  +  »  '  2  ♦  *  T  -  ~'-vr*fr=< 

-i  *  *t  ,,'=*  Flc.c  > 


k  X  *  *  *  y  -  *  ‘ 


. rn«PUTr  rTcFror-.Cr  prTs<rE*‘ 

vqft.'.L^svc  f4.LL4)  '  'J  ’* y 

>  7  '  T;t,4 

'rl  "'I.'"  ■■ !  > 

-fi)=y'<  >-y^(  i*v) 

-  T  r 


’  i  .  L  L*  '  =  v  ~  f  4  •  L  L  4  )  * -J  *  v  ,  , 

__  _  .  _  .  .rv.ruTr  ,<  -r-rtrTrvrc  F.^  I  '•  T  F  ’  L  A  T I 


1  3  T  '1  <  '* 

•f -a  (  t  , :  > 
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